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SUMMARY 


An  investigation  of  the  combustion  properties  of  carbon  agglom¬ 
erates,  formed  during  the  combustion  of  carbon-black  slurry  fuels,  is 
described.  The  carbon  agglomerates  of  several  noncatalyzed  carbon-black 
formulations  were  considered.  This  Included  monodisperse  blacks  having 
ultimate  carbon  particle  sizes  of  70,  150  and  300  Qm  as  well  as  a  blend 
of  the  70  and  300  nm  blacks  (50%  each  by  mass) .  Since  carbon-black 
slurry  fuel  formulation  efforts  are  continuing,  the  formulations  tested 
are  not  necessarily  optimum. 

Observations  of  the  combustion  properties  of  carbon  agglomerates 
were  made  by  injecting  dried  agglomerates  (initial  diameters  of  10-75  urn) 
into  the  post-flame  region  of  a  flat-flame  burner.  The  burner  was 
operated  at  atmospheric  pressure  with  fuel  equivalence  ratios  in  the 
range  0. 2-1.4  and  flame  temperatures  in  the  range  1620-1960  K.  Gas 
compositions  were  measured  by  isokinetic  sampling  and  analysis  with  a 
gas  chromatograph;  gas  temperatures  were  measured  with  a  coated  fine-wire 
thermocouple — corrected  for  radiation  errors;  gas  velocities  were  mea¬ 
sured  using  laser  Doppler  anemometry.  Particle  diameter,  mass,  tempera¬ 
ture,  velocity  and  residence  time  were  measured  as  a  function  of  position 
above  the  point  of  particle  Injection  using  a  combination  of  sampling  and 
non- intrusive  methods. 

The  measurements  were  interpreted  using  an  analysis  of  agglomerate 
combustion  developed  during  earlier  work  with  carbon-black  slurries  in 
this  laboratory.  The  model  provides  predictions  of  particle-life  history 
i.e.,  the  variation  of  diameter,  mass  and  temperature  as  a  function  of 


SUMMARY  (Cont'd.) 


time  in  Che  burner  gases.  The  analysis  treats  heat  and  mass  transfer 
processes  in  the  convective  environment  surrounding  the  carbon  agglom¬ 
erate  as  well  as  radiative  heat  transfer  between  the  agglomerate  and 
its  surroundings.  Two  carbon  reaction  mechanisms  are  considered:  one 
assuming  reaction  between  carbon  and  0^,  CO^  and  ^0;  and  one  assuming 
reaction  between  carbon  and  OH  (high  equivalence  ratios)  or  0 2  (low 
equivalence  ratios).  The  carbon  agglomerates  have  a  porous  structure. 
This  was  treated  using  a  shrinking-sphere  model  but  employing  empirical 
transport  enhancement  and  area/reactivity  multiplication  factors  to 
allow  for  effects  of  pores. 

The  model  provided  a  good  correlation  between  predictions  and  mea¬ 
surements  using  either  reaction  mechanism.  The  empirical  parameters  were 
nearly  independent  of  flame  conditions  and  initial  agglomerate  diameter. 
After  an  initial  period  of  pore  development,  the  empirical  parameters 
were  also  relatively  Independent  of  mass  fraction  of  carbon  that  was 
reacted.  In  general,  the  empirical  parameters  provide  a  useful  method  of 
summarizing  the  properties  of  a  carbon-black  formulation,  since  their 
values  tend  to  increase  for  faster  reacting  materials. 

The  rate  of  combustion  of  monodisperse  blacks  increased  slightly 
as  the  ultimate  carbon-black  particle  size  was  reduced,  i.e.,  size 
reductions  from  300  to  70  nm  yielded  10-20%  reductions  in  particle 
residence  times  to  react  fixed  fractions  of  the  carbon.  The  carbon-black 
blend  had  an  initial  density  that  was  16%  larger  than  the  monodisperse 
agglomerates.  Times  to  react  fixed  fractions  of  the  carbon  in  the  blend 
were  also  longer  than  for  monodisperse  agglomerates;  ranging  from  10-20% 
longer  at  diffusion-controlled  conditions  to  as  much  as  50%  longer  at 
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kinetic-controlled  conditions.  These  increases  are  only  partly  explained 

by  the  increased  density,  suggesting  that  the  presence  of  the  smaller 

particles  in  the  blend  tends  to  inhibit  the  development  of  pores. 

Variations  in  initial  agglomerate  diameter,  d  ,  and  flame  con- 

po 

ditions  influenced  residence  times  required  for  reaction  to  a  much 

greater  degree  than  the  carbon-black  formulation.  Reaction  residence 

times  were  proportional  to  d  ^  t0  for  the  low  particle  Reynolds 

po 

numbers  (0-0.35)  of  the  present  tests — tending  toward  the  lower  power 
near  kinetic-controlled  conditions  and  the  higher  power  near  diffusion- 
controlled  conditions.  Residence  time  requirements  decreased  mono- 
tonically  as  flame  equivalence  ratio  was  reduced,  for  the  relatively 
high  gas  temperatures  of  the  present  tests.  Therefore,  effects  of 
atomization  and  mixing  in  the  combustion  chamber  are  very  significant 
for  achieving  good  combustion  efficiencies  with  carbon-black  alurry  fuels. 


xv 


1.  Introduction 


1.1  General  Objective* 

Liquid  fuels  having  high  energy  densities  (high  energy  release 
when  burned  with  air  per  unit  volune  of  fuel)  are  needed  to  lap rove  the 
perforaance  of  volume- limited  propulsion  systems  [1] .*  It  has  become  in¬ 
creasingly  difficult,  however,  to  obtain  additional  iaproveaents  employing 
liquids  alone  [2].  Therefore,  slurries  are  increasingly  receiving 
attention  since  they  can  be  handled  and  burned  similar  to  liquids,  but 
have  higher  energy  densities  due  to  the  presence  of  solids  [2-5]. 

The  combustion  properties  of  slurry  fuels  are  not  well  known.  Past 
work,  however,  has  shown  chat  combustion  of  the  solid  in  a  slurry  requires 
substantially  longer  combustion  chamber  residence  times  than  conventional 
liquid  fuels  [2-7].  This  implies  that  effective  use  of  slurries  will 
require  modification  of  combustor  designs  from  those  used  for  conventional 
liquid  hydrocarbons.  Slurry  fuel  formulation  also  has  a  significant 
effect  on  residence  time  requirements,  since  solids  loading  and  solid 
particle  properties  affect  both  atomization  and  combustion  properties. 

This  adds  new  dimensions  to  fuel  development  efforts  which  are  not  en¬ 
countered  for  conventional  liquid  fuels. 

The  objective  of  the  present  Investigation  was  to  examine  the  com¬ 
bustion  properties  of  a  particular  class  of  slurry  fuels  of  interest  to 
the  Air  Force,  consisting  of  slurries  using  carbon-black  particles  as  the 
solid  phase.  These  slurries  were  formulated  by  Suntech,  Inc.,  Marcus 
Hook,  PA,  as  part  of  an  Air  Force  sponsored  fuel  development  effort. 

While  the  fuels  investigated  are  representative  of  candidate  carbon  slurry 
fuels,  they  are  not  necessarily  the  most  attractive  formulations  since 

^Numbers  in  brackets  denote  references. 


fuel  development  efforts  ere  continuing.  However,  the  materials  still 
provide  a  reasonable  indication  of  carbon-black  slurry  combustion  prop¬ 
erties,  useful  for  future  fuel  development  efforts. 

The  main  emphasis  of  the  present  investigation  was  to  observe  the 
combustion  of  the  solid  phase  originating  from  individual  slurry  drops. 

In  order  to  simulate  practical  combustor  conditions,  observations  were 
made  of  free  particles  Injected  into  the  post-flame  region  of  a  flat- 
flame  burner  operating  at  atmospheric  pressure  (Initial  particle 
diameters  less  than  100  microns) .  Analysis  of  the  process  was  also  under¬ 
taken  in  order  to  assist  interpretation  of  the  measurements. 

In  the  following,  past  work  on  combustion  properties  of  carbon-black 
slurries  is  reviewed  prior  to  specifically  identifying  the  objectives  of 
the  present  investigation.  Experimental  and  theoretical  methods  used 
during  the  investigation  are  then  described.  The  report  concludes  with 
a  discussion  of  experimental  results  and  the  comparison  between  pre¬ 
dictions  and  measurements.  The  entire  investigation  is  considered  here, 
although  portions  of  the  study  have  also  been  reported  elsewhere  [8-10] . 

1.2  Background 

Existing  information  on  the  combustion  properties  of  carbon- 
black  slurries  is  relatively  limited.  Carbon-black  slurries  have  been 
burned  in  a  gas-turbine  combustor,  by  Bruce  et  al.,  [2,3]  and  in  well- 
stirred  reactors  by  Salveson  [5],  and  Lavld  and  Ruth  [11].  These  studies 
indicated  that  carbon  slurries  required  a  greater  combustor  residence 
time  for  good  combuatioa-eff icieucy- than  conventional  liquid  fuels,  how¬ 
ever,  the  mechanism  of  combustion  was  not  established. 
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Law  and  coworkers  [12,13]  considered  the  related  problem  of  coal 
slurry  combustion — observing  both  supported  and  freely  falling  drops 
burning  In  air.  It  was  found  that  the  process  could  be  broadly  divided 
Into  two  stages.  The  liquid  fuel  evaporated  or  burned  in  the  first 
stage,  leaving  an  irregular  agglomerate  of  the  coal  particles.  The 
agglomerate  subsequently  reacted  or  was  quenched,  depending  upon 
ambient  conditions,  during  the  second  stage  of  the  process. 

Earlier  studies  in  this  laboratory  considered  the  combustion  of 
large  carbon-black  slurry  drops,  400-1000  um  diameter,  supported  at 
various  positions  in  a  turbulent  diffusion  flame  [4,6,7].  It  was 
found  that  carbon-black  slurries  burned  in  a  two-stage  process  similar 
to  coal  slurries.  The  liquid  evaporated  rapidly,  leaving  a  porous  agglom¬ 
erate  containing  all  the  carbon-black  particles  originally  in  the  drop. 

The  agglomerate  then  heated-up  and  either  reacted  or  was  quenched  de¬ 
pending  upon  position  in  the  flame.  Even  at  maximum  agglomerate  reaction 
rates,  however,  the  second  stage  was  at  least  an  order  of  magnitude 
longer  than  the  liquid  evaporation  stage.  Therefore,  agglomerate  com¬ 
bustion  represents  the  rate-controlling  step  for  good  combustion  efficiency. 

A  model  of  slurry  drop  combustion  was  also  developed  during  earlier 
work  in  this  laboratory  [6,7].  Liquid  gasification  was  treated  following 
past  practice  for  analysis  of  sprays  [14-16],  but  extended  to  allow  for 
the  presence  of  solid  in  the  liquid  phase.  Agglomerate  reaction  was 
modeled  by  extending  the  carbon  particle  combustion  analysis  of  Libby  and 
Blake  [17].  Two  carbon  reaction  mechanisms  were  considered:  (1)  the 
approach  of  Neoh  et  al.  [18]  for  soot  combustion,  where  carbon  reaction 
with  OH  dominates  at  high  fuel  equivalence  ratios  while  reaction  with  0^ 


becomes  significant  at  low  equivalence  ratios;  and  (2)  the  approach  of 


Libby  and  Blake  {17]  for  reaction  of  carbon  with  0^  and  CO^,  but  extended 
to  include  reaction  with  H^O.  The  analysis  of  liquid  evaporation  was 
successfully  evaluated  using  data  from  supported  pure  liquid  drops  in 
the  diffusion  flame  [8].  The  agglomerate  combustion  model,  using  either 
reaction  mechanism,  also  yielded  good  comparison  with  measurements  for  a 
wide  range  of  flame  conditions— af ter  specifying  empirical  constants 
to  allow  for  effects  of  catalyst  and  pores  on  gas-phase  transport  and 
reactive  surface  area.  Combining  these  two  models  yielded  good  pre¬ 
dictions  of  slurry  drop-life  histories  (the  variation  of  drop/particle 
size  and  temperature  with  time)  throughout  both  stages  of  the  process  [7] 
Major  limitations  of  this  work,  however,  were  that  drop  sizes  were  large 
in  comparison  to  drops  in  conventional  combustors,  only  a  single  slurry 
formulation  was  studied,  and  effects  of  pores  on  the  density  of  the 
agglomerates  was  not  studied.  This  raises  significant  questions  con¬ 
cerning  the  generality  of  the  results— particularly  the  empirical  parame¬ 
ters — for  practical  combustor  conditions.  The  effect  of  fuel  formulation 
— particularly  carbon-black  type— on  slurry  combustion  properties  was 
also  not  established. 

1.3  Specific  Objectives 

The  present  investigation  extends  the  study  of  carbon-black 
slurry  combustion — considering  particle  sizes  more  representative  of 
practical  combustors  (initial  diameters  of  10-75  ym).  Since  the  liquid 
gasification  stage  is  conventional  [6,7],  the  study  was  limited  to 
agglomerate  combustion  properties.  Freely  moving  agglomerate  particles 
were  observed  in  the  post-flame  region  of  a  laminar  flat-flame  burner. 


yielding  the  variation  of  agglomerate  temperature,  diameter,  mass  and 
velocity  as  a  function  of  residence  time. 

In  addition  to  effects  of  overall  agglomerate  size,  experiments  were 
also  conducted  to  determine  the  combustion  properties  of  various  carbon- 
blacks  with  ultimate  carbon  particle  sizes  in  the  range  70-300  nm. 

Limited  testing  was  also  completed  on  a  blend  consisting  of  carbon  blacks 
of  different  ultimate  carbon  particle  size,  since  blends  may  be  required 
to  yield  slurries  having  acceptable  rheological  and  stability  properties. 

All  the  measurements  were  interpreted  using  the  agglomerate  combustion 
model  developed  earlier  [6,7]. 

2.  Experimental  Methods 

2.1  Apparatus 

A  sketch  of  the  flat-flame  burner  apparatus  appears  in  Fig.  1. 

The  arrangement  consists  of  a  slurry-droplet  generator,  a  convergent-flow 
section,  a  liquid  evaporation  section,  and  a  flat-flame  burner.  The  mono- 
disperse  drops  produced  in  the  drop  generator  are  dried  in  a  heated  section 
of  tubing  leaving  an  agglomerate  of  solid -particles.  The  carbon  agglom¬ 
erates  then  pass  through  the  post-flame  region  of  a  flat-flame  burner  where 
observations  are  made.  A  sampling  probe  placed  above  the  flat-flame  burner 
collects  the  carbon  agglomerates  for  further  analysis. 

The  flat-flaue  burner  was  fueled  with  a  mixture  of  gaseous  nitrogen, 
oxygen,  carbon  dioxide,  carbon  monoxide,  methane  and  hydrogen.  The  mixing 
pressure  was  monitored  with  a  Heise  absolute-pressure  gauge  having  a 
pressure  range  of  0.0-0.41  MPa.  Gas  flow  rates  were  controlled  with  pres¬ 
sure  regulators  and  needle  valves  and  measured  with  rotameters.  All  rota¬ 


meters  were  calibrated  using  wet-test  maters. 


Gas  mixing  was  aided  by  f lowing  the  gases  through  a  1.5  m  length  of 
10  mm  I.D.  tubing.  The  gas  mixture  then  entered  the  mixing  chamber  through 
an  inlet  near  its  bottom.  The  gas  mixing  chamber  was  constructed  from 
a  51  mm  nominal  diameter  schedule  40  stainless-steel  welded  pipe  and  had  a 
length  of  300  am.  To  provide  a  thoroughly  mixed  fuel  gas  with  a  uniform- 
velocity  profile,  several  layers  of  small  beads  separated  by  void  spaces 
were  provided  in  the  <gas  mixing  chamber— as  indicated  in  Fig.  1.  A 
flat  flame  was  stabilized  near  the  top  surface  of  the  burner.  The  flame 
was  approximately  1-2  am  thick  and  was  completely  uniform  and  stable  at 
all  operating  conditions. 

Equivalence  ratios  between  0.2  and  1.4  were  used  during  the  tests. 

The  temperatures  of  these  flames  could  be  varied  between  1600  K  and  the 
softening  temperature  of  stainless  steel.  By  blending  carbon  monoxide, 
carbon  dioxide,  hydrogen  and  methane,  the  C/H  ratio  could  be  varied  over  a 
very  wide  range.  Finally,  the  bed  of  stainless-steel  balls  at  the  top  of 
the  burner  inhibited  flash-back  and  aided  flame  attachment  so  that  gas 
velocities  could  be  varied  between  0.5  and  3.0  m/s. 

Once  the  flame  was  ignited,  a  53  mm  I.D.  quartz  tube  was  placed  on 
the  burner  surface.  The  tube  surrounded  the  flame  to  prevent  entrained 
air  from  mixing  with  the  products  of  combustion  and  lowering  the  *lame 
temperature.  The  quartz  tube  also  helped  maintain  a  uniform  high  tempera¬ 
ture  of  the  products  of  combustion  far  above  the  burner.  The  quartz  tube 
was  cut  into  lengths  of  50  mm  to  200  mm,  in  50  mm  increments,  to  facilitate 
the  access  of  measuring  probes  at  various  locations  above  the  burner 
surface. 

The  drop  generator  consisted  of  a  TSI  model  3050  Berglund-Liu 
vibra ting-orifice  monodisperse-aerosol  generator,  a  Sage  model  355 


contlnuously-variable  syringe  pump,  end  e  Hewlett-Packard  model  3310B 
function  generator.  Three  orifices  were  used  in  the  Berglund-Liu  drop 
generator:  a  TSZ  model  10336  100  \m  orifice,  a  TSI  model  10355  50  pm 
orifice  and  a  TSI  model  10353  20  im  orifice.  Both  a  10  and  50  ml  reusable 
multiflt  syringe  with  metal  luer  tips  were  used  with  the  syringe  pump. 

A  mixture  of  methane,  oxygen,  carbon  monoxide  and  nitrogen  was  used 
in  the  Berglund-Liu  droplet  generator  to  dilute  and  disperse  the  drops. 
These  gases  were  mixed  to  yield  the  same  equivalence  ratio  and  C/H  ratio  as 
the  fuel  mixture  in  the  flat-flame  burner  assembly  after  evaporation  of 
the  hydrocarbon  carrier  of  the  slurry.  The  flow  rates  of  the  dilution  and 
dispersion  flows  were  metered  and  measured  similar  to  the  main  burner  gases 

The  slurry  had  to  be  thinned  in  order  to  generate  slurry  droplets  with 
the  drop  generator  without  clogging  the  orifice.  This  involved  mixing 
JP-10  with  the  slurry  fuel.  It  was  found  that  the  carbon  particles 
remained  suspended  in  the  mixture  for  several  hours.  Thus,  mixtures  of 
carbon  slurry  and  n-pentane  ranging  from  one  to  five  parts  n-pentane  per 
part  of  volume  of  slurry  were  used  in  the  droplet  generating  system. 

The  syringe  pump  was  adjusted  to  the  minimum  velocity  necessary  to 
establish  a  liquid  jet  emerging  from  the  orifice  in  the  Berglund-Liu 
droplet  generator.  The  frequency  of  the  signal  generator  was  set  at  a 
value  between  1000  and  4000  Hz.  Thus,  between  1000  and  4000  drops  per 
second  were  formed  ranging  in  initial  diameter  between  24  and  270  ym. 

The  50  ml  syringe  allowed  continuous  operation  for  up  to  3.5  hours. 

The  Berglund-Liu  droplet  generator  exhausted  the  slurry  droplets 
into  a  50  mm  length  of  a  115  mm  diameter  plexiglass  tube.  The  plexiglass 
tube  was  connected  to  a  stainless-steel  45°  convergent  cone,  which  was 
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attached  to  the  plexiglass  with  silicon  sealant.  The  cone  was  40  imn  high, 
had  a  115  mm  diameter  opening  at  the  base  and  a  3.4  mm  diameter  opening  at 
its  tip.  The  outlet  of  the  cone  was  welded  to  a  tube  having  an  inner 
diameter  of  3.4  mm  and  a  length  of  610  mm.  The  top  end  of  this  tube  was 
welded  to  a  tube  with  a  1.8  mm  inner  diameter  and  a  length  of  460  mm. 

This  tube  passed  through  a  Conax  thermocouple  seal  at  the  bottom  of  the 
mixing  chamber.  The  top  end  of  the  1.8  mm  diameter  tube  was  flush  with 
the  surface  of  the  flat-flame  burner. 

r 

Nlchrome  heater  wire  and  fiberfrax  ceramic  insulation  were  wrapped 
around  the  outside  of  the  3.4  mm  diameter  section.  Thus,  when  a  particle 
left  the  Berglund-Liu  droplet  generator  it  was  accelerated  through  the 
dohvergent  core  section  of  the  droplet  guide.  Next,  the  droplet  entered 
the  3.4  mm  diameter  section,  which  was  heated  to  600  K.  In  this  region 
the  liquid  fuel  evaporated  leaving  behind  an  agglomerate  of  solid-carbon 
particles.  The  evaporating  liquid  also  served  to  further  accelerate  the 
agglomerate.  When  the  gas  flow  emerged  from  top  of  the  1.8  mm  diameter 
tube  it  had  the  same  equivalence  ratio,  C/H  ratio  and  velocity  as  the 
surrounding  gas  of  the  flat-flame  burner.  The  particle  velocity  was  some¬ 
what  lower  than  the  flow  velocity  due  to  the  gravity  effects  on  the  particle. 
The  initial  size  of  the  carbon  agglomerates  leaving  the  1.8  mm  diameter 
tube  ranged  between  10  and  75  ym.  The  particles  continued  to  flow  up 
through  the  quartz  tube  as  they  reacted. 

Figure  2  is  a  time  exposure  photograph  of  the  carbon  agglomerates 
passing  through  the  flat-flame  burner  and  reacting.  The  quartz  tube  shown 
in  Fig.  1  was  not  in  place  for  this  photograph,  and  the  particles  cool  upon 
exiting  the  hot  potential  core  of  the  flame.  Although  the  particles  do 
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disperse  somewhat,  the  photograph  shows  that  the  particles  travel  in  a 


straight  line  producing  a  streak  approximately  3  mm  wide  at  the  base  for 
an  initial  particle  diameter  of  73.4  ym. 

2.2  Particle  Environment  Measurements 

2.2.1  Velocity  Measurements 

Gas  velocities  were  measured  along  the  axis  of  the  burner 
using  a  laser  Doppler  anemometer  (LDA) .  This  involved  a  5  mW  helium-neon 
laser  operating  in  the  dual-beam  backscatter  mode  (241  mm  focal  length, 
free  aperture  of  60  mm,  50  mm  beam  spacing,  photodetector  aperture  of  0.256 
mm,  with  200  mm  focal  length  detector  focusing  lens) .  This  arrangement 
produced  an  ellipsoidal  probe  volume  0.31  mm  in  diameter  and  3  mm  long. 

The  LDA  output  was  processed  using  a  frequency  tracker.  The  gas  flow  was 
seeded  with  0.3  urn  aluminum  oxide  particles  for  velocity  determination* 
using  a  reverse  cyclone  seeder  [19) . 

2.2.2  Temperature  Measurements 

The  gas  temperatures  along  the  axis  of  the  post-flame 
region  of  the  flat-flame  burner  were  measured  using  a  thermocouple  probe. 
The  probe  was  constructed  from  25  ym  diameter  platinum/platinum-10%  rhodium 
wires  which  were  spot-welded  onto  50  ym  diameter  lead  wires  of  the  same 
material.  The  probe  was  mounted  on  a  traversing  mechanism  which  allowed 
it  to  be  positioned  at  various  locations  along  the  flame  axis.  The  emmis- 
ivity  of  the  thermocouple  was  maintained  at  approximately  0.22  +  0.02  by 
coating  the  thermocouple  with  a  layer  of  silica  roughly  3  ym  thick  [20]. 

The  silica  layer  also  eliminated  catalytic  decomposition  heating  of  the 


The  thermocouple  probe  was  calibrated  with  a  50  pn  chromel/alumel 
thermocouple  in  a  furnace  environment  where  the  radiation  errors  are 
minimized.  Both  thermocouples  agreed  within  6°C,  up  to  a  temperature  of 
1500  K,  which  was  the  maximum  temperature  of  the  furnace.  The  maximum 
radiation  correction  for  this  temperature  probe  was  39°C  over  the  test 
range.  The  reference  junction  was  at  the  ambient  temperature,  shielded 
from  the  flame.  An  integrating  digital  voltmeter  was  used  to  average 
the  signal  over  a  two-minute  period. 

2.2.3  Concentration  Measurements 

The  concentration  measurements  were  made  using  a  gas 
chromatograph  (Varian  model  3720)  incorporating  a  thermal-conductivity 
detector.  The  gas  samples  were  extracted  from  the  flow  at  nearly 
isokinetic  conditions  using  a  water-cooled  probe  constructed  of  stainless 
steel.  The  inlet  of  the  probe  was  0.7  mm  in  diameter.  The  temperature 
of  the  cooling  water  was  maintained  constant  about  338  K  (+  2  K)  to  avoid 
the  condensation  of  fuel  and  water  vapor  in  the  probe.  An  ice  bath  cold 
trap  was  placed  in  the  sampling  line,  however,  in  order  to  condense  water 
prior  to  gas  analysis. 

The  samples  were  analyzed  using  a  carbosieve  column  (Supleco,  Carbo- 
sleve,  type  S,  100/120  mesh,  2.1  m  x  3.2  mm  stainless  steel;  4  min.  hold 
at  35°C;  programming,  25°C  to  175°C  at  25°C/min.)  for  the  analysis  of 
hydrogen,  oxygen,  nitrogen,  carbon  monoxide,  carbon  dioxide  and  methane. 
The  concentration  of  water  vapor  was  obtained  by  a  conservation  of  element 
analysis.  Helium  at  a  flow  rate  of  30  ml/min.  was  used  as  carrier  gas. 

The  output  of  the  hot-wire  thermal-conductivity  detector  was  recorded  with 
a  Varian  strip  chart  recorder. 
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The  Carboaieve  column  was  calibrated  for  hydrogen,  oxygen,  nitrogen, 
carbon  monoxide,  carbon  dioxide  and  methane  using  calibration  gases  manu¬ 
factured  by  the  Scott  Company.  A  0.5  ml  sample  size  was  used  for  all 
composition  measurements. 

2.3  Particle  Measurements 

2.3.1  Velocity  Measurements 

The  particle  velocity  measurements  were  made  with  the  LDA 
system  used  to  measure  the  velocity  of  the  gas  above  the  flat-flame  burner. 
The  particle  spacing  was  also  measured  with  the  LDA  system.  A  Nicolet, 
model  206,  digital  oscilloscope  was  employed  to  measure  the  time  between 
the  pedestal  signals  generated  by  individual  particles.  Since  the  particle 
velocity  was  known,  the  distance  between  particles  could  be  determined. 

2.3.2  Temperature  Measurements 

The  temperature  of  the  carbon  particles  as  they  reacted 
above  the  flat-flame  burner  assembly  was  measured  using  an  optical  pyrometer. 
The  pyrometer  used  was  a  Leeds  and  Northrup,  model  8622,  with  a  temperature 
range  of  1000-3150  K.  The  optical  pyrometer  was  calibrated  using  the 
platinum/platinum-10%  rhodium  thermocouple  temperature  probe  employed  for 
gas  temperature  measurements. 

The  particle  temperature  measurement  is  a  brightness  temperature, 
assuming  that  the  particles  are  black  bodies.  This  technique  is  particu¬ 
larly  adapted  for  temperature  measurements  of  high  temperature  particle 
streams,  since  it  is  independent  of  the  number  of  particles  observed  at 
any  instant. 

2.3.3  Diameter  Measurements 

The  probe  used  to  collect  the  carbon-agglomerate  particles 
for  analysis  is  shown  in  Fig.  3.  The  probe  was  positioned  along  the  flame 
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AND  BUBBLE  METER 
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Figure  3.  Sketch  of  the  quenching  probe  used  for  particle  collection. 


centerline  at  various  locations  above  the  flat-flame  burner  with  a  traversing 
mechanism.  The  collection  filter  was  constructed  from  a  Gelman,  model  2220, 
stainless-steel  inline  filter  holder.  The  upstream  end  of  the  filter  holder 
was  modified  to  include  two  concentric  stainless-steel  cones  which  were 
silver  soldered  to  the  upstream  face  of  the  filter  holder.  The  inner  cone 
had  a  6.4  mm  inlet,  a  15°  angle,  and  a  length  of  23.7  mm.  The  outer  cone 
had  a  12.7  mm  inlet,  a  30°  angle,  and  a  length  of  30.5  mm.  Both  cones  had 
an  approximate  thickness  of  1.6  mm.  A  stainless-steel  tip  was  silver 
soldered  to  the  outer  core.  The  tip  was  12.7  mm  in  diameter,  6.4  mm  in 
length  and  tapered  to  a  4.7  mm  inlet  with  an  angle  of  45°. 

A  nitrogen-quenching  flow  passed  through  the  outer  cone  of  the  filter 
holder  from  a  6.4  mm  diameter  tube.  The  nitrogen  passed  between  the  two 
concentric  cones  to  the  inlet  of  the  collection  probe  where  the  reacting 
particles  were  quenched  as  they  entered  the  probe.  The  particles  then  were 
collected  on  47  mm  diameter  Alpha  metricel,  grade  Alpha-200,  Gelman  membrane 
filters,  model  60585.  These  filters  have  a  pore  size  of  0.2  pm,  a  thickness 
of  65  pm,  a  very  high  chemical  resistance,  and  a  maximum  operating  tempera¬ 
ture  of  450  K. 

The  particle  collection  probe  was  connected  to  a  vacuum  pump.  The 
velocity  at  the  inlet  to  the  collection  probe  was  matched  with  the  gas 
velocity  at  the  sampling  location,  using  a  regulating  valve  and  bubble  flow 
meter.  A  sufficient  amount  of  nitrogen  gas  was  mixed  in  the  collection 
filter  to  lower  the  gas  temperature  to  400  K  and  quench  the  reacting 
carbon  particles.  The  operating  pressure  in  the  filter  was  approximately 
71  kPa  which  also  tended  to  help  quench  the  reacting  carbon  particles. 


A  quick  response  toggle  operated  bellows  valve  was  also  placed  in 
the  flow  line  between  the  collection  probe  and  the  vacuum  pump  so  that 
the  collection  filter  could  be  operated  for  short  time  periods.  Small 
sampling  times  resulted  in  widely  spaced  single  particles  on  the  collection 
filter.  After  the  sample  was  collected  the  filter  medium  was  removed  and 
the  sizes  of  the  particles  collected  were  analyzed  on  a  scanning  electron 
microscope  (SEM) . 

The  particle  diameter  quoted  here  is  the  arithmetic  average  of  six 
diameter  measurements  per  particle  each  spaced  30°  apart  and  passing 
through  the  center  of  the  particle.  At  least  eight  particles  were 
examined  in  order  to  obtain  an  average  particle  size. 

2.3.4  Mass  Measurements 

The  filter  collection  probe  illustrated  in  Fig.  3  was 
employed  to  collect  particles  for  particle  mass  measurements.  However, 
the  sampling  period  for  the  mass  measurements  was  increased  to  between 
10  seconds  and  2  hours.  A  stainless-steel  filter  medium  was  also  substi¬ 
tuted  for  the  Alpha  metricel  filter  described  earlier.  The  stainless-steel 
filter  was  a  Dynalloy  filter,  model  X4,  manufactured  by  the  Fluid  Dynamics 
Branch  of  Brunswick  Technetlcs.  The  filters  had  a  diameter  of  47  am,  a 
thickness  of  300  ym,  a  pore  size  of  3  ym.  The  filters  were  preconditioned 
by  baking  them  for  three  hours  at  a  temperature  in  excess  of  1000  K  to  burn 
off  any  contaminants. 

After  the  carbon  particles  were  collected,  the  filter  was  removed  and 
placed  in  an  evacuated  reaction  chamber.  The  temperature  was  raised  to 
500  K,  and  the  pressure  was  lowered  below  4  kPa.  These  conditions  were 
maintained  for  several  hours  to  remove  any  condensibles  entrapped  within  the 
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particles.  The  dried  carbon  particles  were  then  exposed  to  an  oxygen  and 
hydrogen  environment .  The  temperature  was  then  raised  above  1000  K  and 
the  pressure  Increased  to  85  kPa.  These  conditions  were  maintained  for 
several  hours  until  all  of  the  carbon  was  reacted.  By  measuring  the  carbon 
dioxide  and  carbon  monoxide  concentrations  In  the  reaction  chamber,  the 
mass  of  carbon  Initially  present  on  the  filter  was  determined.  Knowing  the 
total  number  of  particles  cbllected  during  the  collection  period  (from  the 
droplet  production  frequency)  gave  the  mass  of  carbon  per  particle.  The 
pressure  of  the  reaction  chamber  was  kept  below  atmospheric  to  monitor 
whether  the  system  developed  leaks,  since  nitrogen  would  show  up  in  the  gas 
analysis  of  the  reacted  gas  if  a  leak  was  present.  The  gas  chromatograph 
system  used  for  these  measurements  was  the  same  as  for  the  particle 
environment  measurements. 

The  carbon  slurry  mixture  incorporated  several  additives  to  assure 
stability.  In  order  to  determine  whether  these  additives  altered  the 
carbon  mass  measurement,  a  slurry  mixture  without  additives  was  prepared 
in  this  laboratory  and  results  with  this  material  were  compared  with  the 
stabilized  samples.  It  was  found  that  mass  measurements  were  independent 
of  the  presence  or  absence  of  additives. 

2.4  Experimental  Conditions 

Experimental  conditions  for  the  present  investigation  are  summar¬ 
ized  in  Table  1.  The  temperature  and  composition  of  major  species  were 
found  from  the  gas  environment  measurements  described  in  Section  2.2.  The 
concentration  of  OH  is  significant  for  one  of  the  carbon  agglomerate 
reaction  mechanisms.  This  concentration  was  found  using  the  NASA  CEC-72 
chemical  equilibrium  computer  program  [21]. 
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The  temperature  and  composition  of  the  burner  gas  were  essentially 
constant  over  the  measuring  region  (0-200  mm  from  the  burner) — except  for 
the  first  2  mm  above  the  burner  surface  where  the  burner  flame  was  present. 
Gas  velocities,  however,  increased  by  as  much  as  a  factor  of  two,  from  the 
values  at  the  burner  surface  given  in  Table  1,  due  to  boundary  layer  de¬ 
velopment  on  the  quartz  tube  surrounding  the  post-flame  region.  A  more 
complete  tabulation  of  burner-gas  properties  is  given  in  Appendix  A. 

Initial  particle  velocities  varied,  depending  upon  size,  due  to 
effects  of  gravity.  The  velocity  defect  of  the  agglomerates,  also  varied 
as  reaction  proceeded.  Therefore,  measured  local  values  of  the  velocity 
defect  were  used  to  estimate  the  transport  properties  of  the  flow  around 
the  particle. 

Agglomerate  spacing  varied  during  the  tests,  reaching  a  minimum  of 
eight  particle  diameters.  Results  of  Sanglovannl  and  Kesten  [22],  for 
monodisperse  streams  of  drops,  suggested  that  interparticle  effects  are 
relatively  small  for  the  present  spacings.  Specific  examination  of 
present  conditions  indicated  that  reaction  of  the  particle  stream  varied 
local  reactant  concentrations  generally  less  than  10%;  therefore,  the  test 
conditions  are  essentially  representative  of  results  for  individual 
particles. 

Measurement  accuracy  varied  significantly  over  the  particle  lifetime. 
Estimates  of  average  errors  in  agglomerate  properties  over  the  test  range 
are:  diameter,  10%;  velocity  defect,  10%;  velocity,  5%;  mass,  28%; 
residence  time,  20%;  axial  location,  5%;  and  temperature,  10%.  Burning 
rates,  based  on  particle  diameter,  had  estimated  errors  of  15-30%  for 
particle  diameters  of  2-75  ym.  All  measurements  were  repeatable  within 
these  ranges. 
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Properties  of  the  carbon-blacks  used  during  the  tests  are  also  sum¬ 
marized  in  Table  1.  They  consisted  of  medium  thermal  carbon  blacks  with 
ultimate  particle  sizes  ranging  from  70-300  nm.  The  fuels  were  originally 
mixed  with  hydrocarbon  carriers  and  stabilized.  However,  they  were  diluted 
for  present  use  in  the  drop  generator  as  noted  earlier.  Since  only  dried 
agglomerates  were  considered  during  the  present  study,  and  no  effect  of 
additives  was  found,  carrier  properties  are  irrelevant  and  have  been 
omitted  from  the  table. 

3.  Theoretical  Methods 

3.1  Description  of  Model 

Slurry  drop  combustion  is  analyzed  In  two  stages.  The  first 
stage  treats  evaporation  of  the  liquid  carrier  to  produce  the  carbon 
agglomerate.  This  portion  of  the  analysis  is  a  straightforward  appli¬ 
cation  of  drop  evaporation  theory,  aside  from  allowing  for  the  thermal 
capacity  and  volume  of  the  dispersed  carbon  particles.  Since  the  liquid 
phase  was  not  considered  during  this  Investigation,  this  analysis  and  its 
evaluation  will  not  be  considered  here.  A  complete  description  of  the 
theoretical  treatment  of  the  liquid  evaporation  period  is  presented  in 
Refs.  6,  8  and  10. 

Given  initial  conditions  based  on  analysis  of  the  liquid  evaporation 
process,  the  agglomerate  combustion  model  provides  the  subsequent  life- 
history  of  the  particle  during  the  second  stage  of  the  process.  This  in¬ 
volves  computation  of  agglomerate  mass,  size,  temperature,  velocity  and 
position  as  a  function  of  time — for  given  ambient  conditions  of  the 
particle.  Only  the  major  features  of  the  model  will  be  presented  here, 
in  order  to  assist  understanding  of  the  present  results.  Complete  details 


of  the  agglomerate  combustion  model  may  be  found  in  Refs.  6,  7  and  10. 

20 


The  model  must  account  for  the  physical  features  of  agglomerates 
formed  from  carbon-black  slurries.  The  nature  of  the  agglomerates  can  be 
seen  from  the  SEM  photographs  of  Fig.  4 — one  at  the  start  of  reaction  and 
a  second  after  12.22  of  the  agglomerate  mass  was  reacted.  The  nonreacted 
particle  has  some  Irregularities,  but  Is  roughly  spherical  with  a  rela¬ 
tively  smooth  surface.  The  particle  diameter  decreases  as  reaction 
proceeds,  however,  the  particle  also  becomes  porous  with  Increased 
porosity  near  the  surface — similar  to  observations  of  large  agglomerates 
[6,7]. 

The  variation  of  particle  density  with  extent  of  reaction  provides  a 

more  quantitative  indication  of  the  extent  of  subsurface  reaction  and 

porosity.  This  is  illustrated  in  Fig.  5  for  the  complete  range  of  particle 

environments  considered  here  (300  nm  ultimate  carbon  particle  size) .  The  ua 

3 

reacted  density  of  1820  kg/m  agrees  with  the  value  reported  by  Bruce  et  al. 
[2]  for  dried  slurries  of  the  same  carbon  black.  The  variation  of  density 
with  extent  of  reaction  is  relatively  independent  of  particle  diameter  and 
flame  condition,  yielding  the  following  empirical  correlation 

P„  »  Pnn  (1  -  E)0*6  (1) 

P  PO 

A  characteristic  thickness  (8)  of  the  pore  zone  can  be  defined  as  the 
difference  in  radius  between  the  particle  and  a  particle  having  the  same 
mass  but  the  initial  density.  The  density  correlation  of  Eq.  (1)  yields 

5/rp  -  1  -  (1  -  £)°* 2  (2) 
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Figure  4.  SEM  photographs  of  collected  particles:  <J>  «  1,  Tf  =  1683  K, 

dpo  “  73.4  ym.  Upper  and  lower  photographs  havereacted  mass 
fractions  of  0.000  and  0.122,  respectively,  and  d  =  300  nm. 
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Figure  5. 


Apparent  particle  density  as  a  function  of 
reacted  mass  fraction  for  d  "  300  nm. 


Equation  (2)  indicates  that  <5/rp  <  25%  for  e  <  75%.  Therefore,  the  particle 
decreases  in  size,  accompanied  by  the  development  of  a  pore  structure  which 
is  confined  to  a  relatively  small  fraction  of  the  current  particle  radius, 
over  much  of  the  particle  lifetime. 

Pore  combustion  models  have  been  constructed  for  burning  coal  particles, 
however,  such  models  are  relatively  complex.  A  great  deal  of  information 
on  pore  structure  is  also  needed  prior  to  computations  with  pore  models. 
While  such  developments  are  a  desirable  long-term  objective,  work  along 
these  lines  was  deferred  for  the  present,  in  favor  of  a  semi-empirical 
treatment  of  effects  of  pore  structure.  A  shrinking  sphere  model  was  con¬ 
sidered,  since  this  approach  represents  the  major  features  of  the  process. 
Empirical  parameters  were  defined  to  treat  enhancement  of  convective 
transport  rates  and  Increased  surface  area  for  reaction  due  to  pores. 

The  latter  parameters  also  represent  effects  of  catalyst  [6,7].  Fortui¬ 
tously,  the  empirical  parameters  are  relatively  independent  of  particle 
size  and  flame  environment,  similar  to  density.  Therefore,  they  provide  a 
convenient  method  for  summarizing  effects  of  pores. 

The  model  was  an  extension  of  the  carbon  combustion  analysis  of  Libby 
and  Blake  [17].  Major  assumptions  of  the  model  are:  spherical  isolated 
particle;  uniform  particle  properties  at  each  instant  consisting  of  pure 
carbon;  quasisteady  gas  phase;  mass  diffusion  only  by  concentration 
gradien's;  all  species  assumed  to  have  equal  molecular  weights  and  specific 
heats;  and  constant  average  gas-phase  properties  at  each  instant  of  time 
(evaluated  at  the  reference  state  given  elsewhere  [6,7]),  and  radiation  from 
the  agglomerate  to  the  surroundings  treated  assuming  a  transparent  flame  and 
an  opaque  blackbody  agglomerate. 
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3.2  Gas-Phaae  Transport 


Since  binary  diffusivities  are  equal,  mass  transport  is  con¬ 
veniently  expressed  in  terms  of  element  mass  fractions.  Four  elements, 
C,  0,  N  and  H,  and  seven  species,  02»  N,,  CO^,  CO,  H^O  and  OH,  appear 
in  the  analysis.  The  element  and  species  mass  fractions  are  related  as 
follows : 


7 

Y  -  £  u..  Y  ,  i  -  1,4  (3) 

j-1  3  3 

where  is  the  mass  fraction  of  the  i'th  element  in  the  j'th  species. 
Given  the  Y^  at  a  particular  location,  Eq.  (3)  provides  four  equations 
to  determine  the  seven  unknown  Y^.  Three  additional  equations  were 
obtained  from  the  assumption  of  thermodynamic  equilibrium.  The  equilib¬ 
rium  relations  used  were: 


K 


Pi 


2H20 


2H_  +  0, 


2H20 


+  20H 


(4) 


2C02  2C0  +  02 

where  the  are  known  functions  of  temperature. 

The  Y^  are  found  by  solving  gas-phase  transport  equations.  During  this 
solution,  effects  of  convection  were  evaluated  using  film  theory,  where  the 
gas  phase  is  treated  as  a  spherically  symmetric  layer  having  an  outer  radius 
found  from  empirical  convection  correlations  for  spheres  as  follows  [23]: 


Vrp  “  Nu/<Nu"2)  5  Vrp  *  Sh/  (Sh-2) 


(5) 
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The  effect  of  flow  through  the  pore  structure  was  'considered  by  defining  a 
transport  enhancement  factor,  where 


Nu  -  ?  Nu’  ;  Sh  -  5  Sh’  (6) 

Nu'  and  Sh'  are  the  Nusselt  and  Sherwood  numbers  of  smooth  spheres  found 
similar  to  past  work  on  drop  transport  properties  [14-16,23]. 

Denoting  carbon  as  element  1,  the  film  theory  analysis  yields  the 
following  expressions  for  the  at  the  particle  surface: 

Yx  *  1  -  (1  -  Y  )  exp  (-2K/Sh)  (7) 

p  Aoo 


Given  the  ambient  conditions,  the  temperature  of  the  particle,  and  K 

Eqs.  (3)-(10)  provide  the  Y^  ,  Y^  and  q”«  The  complete  solution  requires 

P  P 

a  relationship  between  K  and  the  Y^.  This  is  provided  by  the  surface 
reaction  properties  discussed  next. 

3.3  Surface  Reaction  Models 

There  are  numerous  uncertainties  in  modeling  the  reaction  of  a 
porous  carbon  agglomerate  in  a  flame.  The  flame  environment  contains 


several  species  which  are  potential  oxidants  of  carbon,  e.g.,  H^O,  CO^,  0, 
0,f  and  OH,  and  the  reaction  mechanism  Is  not  well  established.  Since  the 
material  Is  porous,  reaction  Is  not  limited  to  the  apparent  outer  surface 
of  the  agglomerate.  Carbon  substances  exhibit  Intrinsic  variations  in 
reactivity  due  to  their  surface  structure.  Finally,  the  presence  of  a 
catalyst  is  known  to  influence  reaction  rates  [6,7]. 

The  effect  of  various  gaseous  reactants  was  treated  by  considering 
a  mechanism  recently  proposed  by  Neoh  et  al.  [18],  as  well  as  an  ex¬ 
tension  of  an  earlier  approach  employed  by  Libby  and  Blake  [17].  The 
effect  of  pores,  surface  reactivity  and  catalysts  was  treated  by  identify¬ 
ing  empirical  area/reactivity  factors,  selected  to  best  match  the  present 
measurements.  (It  was  found  that  a  fixed  value,  for  a  given  carbon  black 
and  reaction,  could  correlate  the  measurements  over  the  range  of  the  data.) 

Neoh  et  al.  [18]  concluded  that  OH  is  the  dominant  carbon  oxidant 
under  fuel-rich  conditions;  that  0  is  of  secondary  importance  at  tempera¬ 
tures  below  2000  K;  and  that  O2  becomes  a  significant  factor  for  fuel-lean 
conditions,  particularly  at  lower  temperatures.  Present  test  conditions 
correspond  to  these  circumstances,  therefore,  carbon  reaction  with  OH  and 
0^  was  considered,  while  ignoring  the  remaining  species.  Neoh  et  al.  [18] 
found  that  reaction  of  carbon  with  OH  could  be  represented  by  assuming  a 
constant  collision  efficiency,  yielding  the  following  reaction  rate 
expression 


Krl  P0H 


(11) 
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The  results  of  Nagle  and  Strickland-Constable  [24]  were  employed  to  de¬ 
termine  the  rate  of  reaction  of  carbon  with  02*  These  rate  expressions 
are  also  in  reasonable  agreement  with  later  measurements  by  Park  and 
Appleton  [25].  The  carbon  reaction  rate  is  given  by 

*2  '  Kr2  P02  X/U  *  Kr3  Po2>  +  Kr*  P02  <X  '  X>  (12) 

where 

X  =»  (1  +  Kr5/(Kr4  P^)]'1  (13) 

The  measurements  of  Neoh  et  al.  [18]  suggest  that  the  larger  of  or 
should  represent  the  reaction  rate  at  any  condition.  Therefore,  the 
dimensionless  mass  burning  rate  was  determined  from 


K  -  (rp/pD)  max  (a.^,  a2R2^ 


(14) 


where  the  a^  appearing  in  Eq.  (14)  are  the  empirical  area/reactivity 
multiplication  factors. 

The  carbon  reaction  mechanism  used  by  Libby  and  Blake  [17]  considers 


reaction  with  0^  and  CO^  employing  the  following  expressions: 


R„  *  K  ,  ?n 
3  r6  0, 


R4  *  Kr7  PC0„ 


(15) 

(16) 


This  approach  was  extended  during  the  present  study  to  include  the  re¬ 
action  of  carbon  with  H^O,  employing  the  results  of  Johnstone  et  al. 


[26].  The  reaction  rate  expression  in  this  case  is 


R5  “  Kr8  PH20/(1  +  Kr9  PH2  +  KrlO 


(17) 
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For  this  second  approach,  Che  total  reaction  rate  of  carbon  was  obtained 
as  the  sun  of  the  rates  for  0^,  C0?  and  H^O,  ignoring  potential  inter¬ 
actions  between  reactants.  This  yields: 

5 

K  -  (r  /pD)  l  a.  R.  <18> 

P  i"3  1  1 

The  specific  reaction  rate  parameters  used  in  the  computations  are 
summarized  in  Table  2,  where  the  are  assumed  to  have  the  following 
general  form: 

ni 

Kri  -  T  1  exp(-E/RT)  (19) 

In  all  cases,  was  found  from  the  mass  fractions  at  the  particle  surface 
as  follows: 

Pi  -  M  Yt  p/Mt  (20) 

The  a^  and  5  were  selected  to  fit  the  data  as  well  as  possible.  Par¬ 
ticular  values  of  these  parameters  will  be  presented  later.  Given  the 

gas-phase  solution  for  a  particular  value  of  K,  the  Y.  are  found  as 

ip 

described  earlier.  Equations  (11)-(14)  for  the  Neoh  et  al.  [18]  approxi¬ 
mation,  or  Eqs.  (15)-(17)  for  the  extended  Libby  and  Blake  [17]  approxi¬ 
mation  yield  a  prediction  of  K,  which  could  then  be  iterated  to  find 
particle  transport  rates. 

3.4  Particle-Life  Histories 

Given  the  temperature  and  size  of  the  particle,  and  the  ambient 

conditions,  the  gas  phase  and  surface  reaction  analyses  yield  the  Y  , 

P 

Y.  ,  K  and  q".  The  variation  of  the  size  and  temperature  of  the  particle 

*p  c 

was  determined  by  solving  equations  for  conservation  of  particle  mass  and 


energy . 


Table  2.  Susnary  of  Reaction  Rate  Parameters  In  K 


rl 


1 

A 

n 

E  (kcal/gmol) 

la 

361  kg  K1/2/m2s  atm 

-1/2  0 

2 

2400  kg/m2 a  atm 

0 

30.0 

3 

21.3  atm"1 

0 

-  4.1 

4 

2 

0.535  kg/m  a  atm 

0 

15.2 

5 

18.1  x  106  kg/m2s 

0 

97.0 

6 

2 

87100  kg/m  s  atm 

0 

35.8 

7 

2 

2470  kg/m  s  atm 

0 

41.9 

8b 

15.15  x  10"3  kg/m2a 

atm  0 

32.7 

9 

9.42  x  10"11  atm"1 

0 

-60.8 

10 

-16  -1 
7.07  x  10  10  atm  A 

0 

-79.3 

aAssumlng 

a  collision  efficiency  of  0.28. 

b Assuming  a  surface  area  of  reaction  of  1.15 

2 

m  /g  of  carbon,  which  is 

the  average  surface  area  over 

the  period  of 

reaction. 

Neglecting  the  relatively  slow  density  variation  of  the  particle,  con¬ 


servation  of  particle  mass  yields: 


d  r 


-pDK/p  r 
P  P 


where  p^  is  found  from  Eq.  (1) 

Conservation  of  energy  at  the  particle  surface  yields: 


(21) 


d  T 


-  i  )  +  q"  + 
c  c 


(22) 


where  q^  is  the  particle  surface  heat  flux  for  radiation  to  the  enclosure 
surrounding  the  flow: 


4  4 

e  (T  -  T  ) 
P  P  w 


(23) 


The  initial  conditions  for  Eqs.  (21)  and  (22)  were  based  on  conditions 
at  the  end  of  the  particle  injection  tube.  At  this  point,  all  the  carbon 
initially  in  the  slurry  was  present  in  the  agglomerate  and  the  agglomerate 
density  and  temperature  are  known. 

Equations  (21)  and  (22)  were  numerically  integrated  on  the  computer, 

using  Gear's  method.  The  nonlinear  algebraic  set  of  equations  which  must 

be  satisfied  in  order  to  determine  the  Y.  ,  Y,  ,  q"  and  K,  were  solved 

lie 
P  P 

using  the  Newton-Raphson  method,  given  the  current  values  of  T^,  r^  and 
the  ambient  conditions. 

Measured  particle  and  gas  velocities  were  used  during  these  compu¬ 


tations.  Residence  time,  position  and  velocity  of  the  particles  were 
related  by  integrating 


Results  and  Discussion 


4.1  Agglomerates  of  Monodisperse  Carbon  Blacks 

4.1.1  Appearance 

The  appearance  of  the  carbon  agglomerates  varied  with 
extent  of  reaction,  but  was  relatively  independent  of  flame  conditions  or 
ultimate  carbon  particle  size  in  the  slurry.  Therefore,  the  SEM  photo¬ 
graphs  in  Fig.  4,  for  du  *  300  nm,  are  representative  of  the  appearance 
of  agglomerates  at  and  near  the  start  of  reaction.  Similar  SEM  photo¬ 
graphs  of  particles  when  the  extent  of  reaction  is  nearly  unity  appear 
in  Fig.  6.  Comparing  all  four  photographs  indicates  that  the  pore 
structure  continues  to  develop  throughout  reaction,  eventually  reaching 
the  center  of  the  particle.  Near  the  end  of  reaction,  the  pore  structure 
develops  to  the  point  where  the  particle  has  a  lacey  structure,  with  an 
apparent  density  near  zero,  but  with  a  substantial  diameter.  This 
behavior  is  reflected  by  the  density  correlation  of  Eq.  (1). 

Similar  to  agglomerate  appearance,  density  was  primarily  a  function 
of  extent  of  reaction.  Results  for  agglomerates  where  d^  »  300  nm  for  a 
wide  range  of  flame  conditions  have  already  been  presented  in  Fig.  5. 
Similar  measurements  for  various  du  are  illustrated  in  Fig.  7,  along 
with  the  correlation  of  Eq.  (1).  These  findings  show  that  Eq.  (1)  provides 
an  adequate  representation  of  agglomerate  density  for  monodisperse  carbon- 
blacks  having  du  in  the  range  70-300  nm. 

4.1.2  Empirical  Parameters 

The  pore  structure  enhances  convective  transport  since  the 
ambient  gas  can  percolate  through  the  porous  particle.  Transport  enhance¬ 
ment  factors  found  from  present  data  for  du  *  300  nm  are  illustrated  in 
Fig.  8.  Results  from  the  large-particle  tests  (400-1000  urn  agglomerates). 


Figure  6.  SEM  photographs  of  collected  particles:  $  ■  1,  *  1683  K 

d  -  73.4  pm.  Upper  and  lower  photographs  have  reacted 
m§ss  fractions  of  0.902  and  0.999  and  d..  “  300  nm. 
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Figure  8. 


Transport-enhancement  factor  as  a  function 
of  reacted  mass  fraction  for  d  *  300  nm. 


for  Che  same  d^,  are  also  shown  on  Che  figure.  This  value  was  taken  Co 
be  conscanc  due  Co  insufficienC  data.  PresenC  tests  considered  particles 
in  the  range  0-75  pm. 

Similar  Co  the  large-particle  tests*  both  energy  and  mass  transport- 
enhancement  factors  were  Identical.  At  the  start  of  reaction,  the 
particles  are  relatively  smooth  and  the  transport-enhancement  factor  is 
nearly  unity.  As  reaction  proceeds,  however,  the  pore  structure  develops 
and  the  transport  enhancement  factor  Increases.  The  final  reduction  in 
£  is  associated  with  the  pore  zone  reaching  the  center  of  the  particle. 

For  £  >  20%,  C  is  within  +  20%  of  the  value  found  earlier  for  large 
particles  [6,7],  where  effects  of  £  were  not  resolved.  This  supports 
the  conclusion  that  the  effect  of  agglomerate  size  on  C  is  not  large. 

Measurements  of  transport  enhancement  factors  for  carbon-blacks  with 

du  in  the  range  70-300  pm  are  illustrated  in  Fig.  9  for  various  flame 

conditions*  These  results  indicate  that  effects  of  d  on  {  are  not 

u 

large  for  this  range  of  ultimate  carbon  particle  sizes. 

Area/reactivity  factors  (a^  for  the  °2“C02_H2°  reaction  mechanism  are 
illustrated  in  Fig.  10  for  dy  -  300  nm.  The  results  are  plotted  as  a 
function  of  reacted  mass  fraction  for  various  flame  conditions  and  initial 
agglomerate  diameters.  The  results  found  for  large  supported  agglomerates 
are  also  shown  on  the  figure  [6,7],  although  effects  of  e  were  not  resolved 
for  this  data.  The  a^  exhibit  less  variation  with  e  than  c,  indicating 
that  microroughness  which  provides  reaction  surface  is  less  influenced  by 
the  development  of  pores.  However,  the  a^  do  increase  somewhat  as  the 
pore  structure  develops.  Present  values  of  the  a^  are  comparable  to  the 
large  particle  values.  The  greatest  discrepancy  between  the  present 

* 

Results  for  dy  ■  300  nm  are  presented  as  a  data  correlation. 
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REACTIVITY  MULTIPLICATION  FACTORS 


measurements  and  Refs.  6  and  7  appears  for  a^  q.  This  difference  is  not 
significant,  however,  since  a„  n  was  arbitrarily  set  equal  to  a„_  during 

“2^  Cu 2 

the  earlier  work,  since  C/H  ratios  were  not  varied  during  these  tests. 

The  present  measurements  for  various  C/H  ratios  indicate  the  a^  q  should 

be  somewhat  greater  than  a_n  . 

cu 2 

The  effect  of  ultimate  carbon  particle  size  on  the  a^  for  the 

O^-CO^-H^O  reaction  mechanism  is  illustrated  in  Fig.  11.  Decreasing  du 

causes  a  progressive  increase  in  the  a^,  for  all  values  of  £.  This  trend 

is  reasonable,  since  smaller  ultimate  carbon  particle  sizes  should  yield 

higher  reactive  surface  areas  per  unit  volume. 

The  a^  for  the  OH-O^  reaction  mechanism  are  Illustrated  in  Fig.  12 

for  d^  =  300  nm.  The  large  particle  values  for  the  same  du  are  also  shown 

on  the  figure.  The  a±  are  also  relatively  independent  of  flame  conditions 

in  the  case,  tending  to  increase  with  increasing  £  and  decreasing  d^Q. 

The  present  values  are  within  20%  of  the  values  found  earlier  for  large 

particles — where  effects  of  e  were  not  resolved. 

The  effect  of  d  on  the  a.  for  the  OH-O.  reaction  mechanism  is  illus- 
u  1  2 

trated  in  Figure  13.  Similar  to  the  02”^2~^20  mechanism,  the  a^  pro¬ 
gressively  increase  as  the  ultimate  carbon  particle  size  decreases. 

Agglomerate  densities,  £  and  the  a^  were  correlated  as  functions  of 
e,  d^and  d^  over  the  present  test  range.  Since  effects  of  flame  conditions 
were  small,  they  were  ignored  in  the  correlations.  These  correlations  are 

summarized  in  Table  3.  Effects  of  d  and  d  ,  while  measurable,  are  not 

p  u 

large.  Except  for  particle  density,  the  effect  of  £  is  also  not  large  for 
£  >  0.20.  Therefore,  mean  values  of  £  and  the  a^  are  also  summarized  in 
Table  3  for  the  present  data,  along  with  mean  values  found  earlier  for 


large  particles.  The  values  of  the  two  sets  of  measurements  are  not  much 
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Figure  12.  Area/reaccivity  multiplication  factors  for  the 
OH-0,,  reaction  mechanism  for  d  ■  300  nm. 
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Szekely  and  Faeth  [6,7],  400-1000  M<n  diameter  supported  particles,  particle  Reynolds  numbers  of 
10-200,  du  -  300  nm.  Ho  variation  of  parameters  was  found  in  this  case  due  to  insufficient  data. 

^Initial  densities  are  1820  (du  -  300  nm),  1859  (d  -  150  nm)  and  1933  (d  -  70  nm)— all  in  kg/m3. 


different,  supporting  the  conclusion  that  the  effect  of  dy  and  d^  on  reaction 
rates  is  relatively  small  for  carbon  black  agglomerates.  The  a^  are  all 
on  the  order  of  100,  which  is  reasonable  in  view  of  the  large  pore  area  per 
unit  surface  area  seen  in  Figs.  4  and  6. 

4.1.3  Burning  Rates 

In  the  following,  we  examine  the  properties  of  agglomerate 
burning  rates  as  a  function  of  flame  condition,  agglomerate  diameter  and 
ultimate  carbon  particle  size.  The  data  are  summarized  in  Appendix  A. 
Predictions  employ  the  empirical  parameters  summarized  in  Table  3. 

Predicted  and  measured  burning  rates  for  agglomerates  having  du  ■  300  nm 
are  illustrated  in  Figs.  14-18.  In  these  figures  particle  burning  rate 
(based  on  rate  of  change  of  the  particle  diameter  with  time)  is  plotted  as 
a  function  of  particle  diameter  with  either  flame  temperature,  equivalence 
ratio  or  initial  particle  diameter  taken  as  a  parameter.  For  these  tests, 
the  burning  rate  based  on  particle  diameter  had  an  estimated  error  of 
less  than  15%  for  a  75  pm  particle  and  less  than  30Z  for  a  2  ym  particle. 
Three  sets  of  predictions  are  shown  using  the  OH-O^  and  02-C02-H20  reaction 
models  with  the  complete  correlations  for  empirical  parameters  given  in 
Table  3,  as  well  as  the  latter  correlation  with  the  constant  average  values 
of  empirical  parameters  given  in  Table  3. 

Predictions  for  either  reaction  mechanism,  using  complete  correlations 
of  empirical  parameters,  are  in  reasonably  good  agreement  with  each  other 
and  with  the  data  in  Figs.  14-18.  Maximum  discrepancies  are  on  the  order 
of  20%.  Predictions  using  constant  average  values  of  the  empirical 
parameters  are  only  slightly  less  satisfactory.  This  is  partly  expected 
since  the  parameters  were  selected  to  fit  the  data,  however,  it  is  most 
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Figure  14.  Variation  of  burning  rate  with  particle  diameter 
and  temperature  for  $  ■  0.6  and  du  ■  300  nm. 
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Figure  18.  Variation  of  burning  rate  with  particle  diameter  and 
equivalence  ratio  for  ■  1950  K  and  d^  *  300  no. 
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encouraging  Chat  Che  analysis  can  treat  effects  of  flame  conditions  and 
particle  size  with  this  single  set  of  parameters. 

Agglomerate  burning  rates  are  independent  of  diameter  under  kinetic- 
controlled  conditions,  where  reaction  rates  are  low  in  comparison  to  rates 
of  mass  diffusion.  For  diffusion-controlled  conditions,  where  reaction 
rates  are  high  in  comparison  to  diffusion  rates,  the  burning  rate  is 
inversely  proportional  to  d  to  where  the  higher  power  is  more 

representative  of  present  tests  at  relatively  low  Reynolds  numbers  [6,7]. 

The  results  at  high  equivalence  ratios,  in  Figs.  15-18,  approach  klnetic- 
controlled  conditions — only  tending  toward  diffusion  control  at  the 
highest  temperatures.  The  low  equivalence  ratio  results  of  Fig.  14 
approach  diffusion  control — except  at  the  lowest  flame  temperatures.  In 
all  cases,  the  process  tends  toward  kinetic-control  for  the  smallest 
agglomerate  diameters,  since  small  particles  provide  little  resistance 
to  convective  heat  and  mass  transfer.  The  very  high  burning  rate  at 
i)i  *  0.2  in  Fig.  16  is  somewhat  misleading,  since  it  is  difficult  to  approach 
this  temperature  level  with  a  conventional  adiabatic  flame,  e.g.,  <t>  ■  0.2 
was  a  quenching  condition  for  the  adiabatic  flame  tests  Refs.  6  and  7. 

Burning  rate  results  for  agglomerates  having  du  -  70  and  150  nm  were 
generally  similar  to  the  findings  pictured  in  Figs.  14-18.  As  dy  is  re¬ 
duced,  C  increases  slightly  while  the  a^  generally  increase.  The  net 
effect  is  a  slight  increase  in  burning  rates  near  diffusion-controlled 
conditions  and  a  slight  increase  in  burning  rates  near  kinetic-controlled 
conditions.  These  trends  are  reasonable,  since  smaller  ultimate  particle 
sizes  are  likely  to  modify  the  effective  roughness  of  the  surface,  which 
causes  variations  in  transport  enhancement,  while  increasing  the  reaction 
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area  per  unit  volume  causing  the  a±  to  increase.  Based  on  this  observation, 
it  seems  likely  that  a  pore  model  could  be  constructed  to  predict  such 
trends  from  the  fundamental  properties  of  the  slurry  formulation. 

4.1.4  Life  Histories 

A  more  complete  picture  of  the  combustion  properties  of 
carbon-black  agglomerates  can  be  seen  from  predicted  and  measured  agglom¬ 
erate-life  histories.  Results  along  these  lines  will  be  considered  in 
some  detail,  since  they  provide  a  stringent  test  of  the  model  as  well  as  a 
survey  of  effects  of  particle  size  and  flame  conditions  on  particle  life¬ 
time.  Naturally,  these  results  are  not  completely  indicative  of  agglom¬ 
erate  behavior  in  an  actual  combustor.  Under  combustor  conditions, 
ambient  properties  of  the  particles  change  continuously,  while  present 
tests  involved  a  fixed  ambient  environment. 

Measured  agglomerate-life  histories  for  du  *  300  nm  are  illustrated 
in  Figs.  19-32  for  various  initial  agglomerate  diameters  and  flame 
conditions.  The  data  involves  particle  diameter  and  mass,  normalized  by 
their  initial  values,  and  particle  temperature,  normalized  by  the  flame 
temperature,  as  a  function  of  time  after  injection  into  the  burner  gas. 
Predictions  for  both  reaction  mechanisms  are  also  shown  on  the  figures, 
using  the  correlations  for  empirical  parameters  given  in  Table  3.  Use  of 
constant  average  empirical  parameters  is  not  considered  here,  since  the 
high  value  of  £,  representative  of  e  >  0.20,  tends  to  overestimate 
particle  heat-up  times. 

The  predictions  of  both  reaction  mechanisms  and  the  measurements  are 
in  reasonably  good  agreement  in  Figs.  19-32.  This  is  partly  expected  since 
the  model  involved  correlation  of  burning  rate  measurements  over  the  same 
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Figure  20.  Particle-life  history  for  $  ■  0.6, 
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Figure  22.  Particle-life  history  for  <j>  *  0.6, 
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Figure  24.  Particle-life  history  for  $  -  1.0, 
T,  ■  1683  K,  d  ■  300  run  and  d  * 
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Figure  27.  Particle-life  history  for  $  -  1.0, 

T.  ■  1819  K,  d  ■  300  rua  and  d _ ■  73.4  pm. 
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range  of  conditions.  However,  the  results  also  indicate  that  the  model 
is  effective  in  predicting  particle  heat-up,  where  reaction  effects  are 
small,  and  final  particle  temperature  levels,  where  radiation  and  chemical 
energy  release  rates  are  important. 

The  time  of  particle  heat-up  is  relatively  short  in  comparison  to 
the  particle  lifetime.  Maximum  particle  temperatures  exceed  local  gas 
temperatures  when  the  reaction  rate  is  high,  e.g.,  $  -  0.6  in  Figs.  20-22; 
and  are  near  or  slightly  below  gas  temperatures,  due  to  radiation  heat 
losses,  when  reaction  rates  are  low,  e.g..  Figs.  19  and  23-32.  Differences 
between  gas  and  particle  temperatures  after  heat-up  are  smaller  than  for 
the  large  particles  observed  earlier  [6,7],  since  the  small  particle  size 
of  the  present  tests  enhances  convective  transport.  Due  to  the  inter¬ 
action  of  diffusion  and  reaction  effects,  agglomerate  lifetime  is  a 
complex  function  of  flame  conditions  and  initial  size.  The  main  generality 

1  t*o  2 

that  can  be  made  is  that  the  lifetime  is  proportional  to  dpQ  ,  where 

the  latter  power  applies  to  diffusion-controlled  conditions.  Therefore, 
good  slurry  atomization  is  required  to  achieve  high  combustion  efficiencies 
within  a  reasonable  combustion  chamber  volume.  Lifetimes  are  shortest  for 
<p  -  0. 6-0.8  (considering  present  findings  as  well  as  Refs.  6  and  7);  there¬ 
fore,  particle  residence  times  at  these  conditions  should  be  maximized  to 
achieve  high  combustion  efficiency  in  a  minimum  combustor  volume. 

Typical  of  any  gasifying  spherical  particle,  normalized  particle  mass 
decreases  much  more  rapidly  than  normalized  diameter;  however,  the  rate  of 
mass  reduction  is  even  greater  for  burning  agglomerates  due  to  the  reduction 
in  particle  density  as  the  pores  develop.  Therefore,  conditions  yielding 
a  relatively  high  mass  burnout  of  carbon  could  still  yield  relatively  large 
agglomerate  particles  at  the  exit  of  the  combustor. 


66 


Particle-life  histories  for  agglomerates  formed  from  carbon-blacks 
having  d^  =*  70  and  150  nm  are  illustrated  in  Figs.  33-38.  Flame  conditions 
for  these  results  were  chosen  to  represent  near  kinetic-controlled  (Figs. 

37  and  38),  near  diffusion-controlled  (Figs.  35-36)  and  intermediate 
(Figs.  33-34)  conditions.  Predictions  of  both  reaction  models  are  also 
plotted  on  the  figures,  using  the  correlations  for  empirical  parameters 
appearing  in  Table  3.  In  general,  the  comparison  between  predictions  and 
measurements  '.s  similar  to  the  results  for  d^  *  300  nm  appearing  in  Figs. 
19-32. 

The  effect  of  ultimate  carbon  particle  size  on  the  agglomerate-life 
history  varied  depending  on  flame  conditions.  Trends  with  respect  to  d^ 
are  most  easily  seen  by  comparing  residence  times  required  to  react  fixed 
fractions  of  the  initial  agglomerate  mass.  These  times  are  summarized  in 
Table  4  for  test  conditions  where  all  three  values  of  dy  were  considered. 
Measured  residence  times  to  react  70  and  90%  of  the  carbon  mass  (t^^  and 
t90%^  are  s*lown  *-n  c^e  table.  Theoretical  values  are  nearly  the  same, 
except  for  <J>  *  0.6,  T^  *  1663  K  and  d^  ■  300  nm,  where  the  measured  values 
are  significantly  larger  than  predicted,  in  which  case  both  values  are 
shown  in  the  table.  Since  the  theory  provides  satisfactory  predictions 
for  the  rest  of  the  data  base,  this  particular  measurement  may  be  anomalous. 
The  plot  of  the  data  appearing  in  Fig.  20  for  this  condition  also  suggests 
problems  with  this  measurement. 

In  general,  reducing  du  reduces  residence  times  required  for  carbon 
consumption.  Aside  from  the  condition  mentioned  above,  when  the  prediction 
is  probably  the  best  indication  of  residence  times,  residence  times  are 
reduced  in  the  range  10-20%  a~  du  is  decreased  from  300  to  7C  nm.  The 
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Figure  33.  Particle-life  history  for  4>  m  0.6,  Tf  *  1663  K,  d  =  150  nm  and  d  »  73.4  pm. 
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Figure  34.  Particle-life  history  for  <P  m  0.6,  Tf  ■  1663  K,  d  ■  70  nm  and  d  *  73.4  pm. 


Particle-life  history  for  ♦  -  1.0,  T-  -  1683  K,  d  -  150  no  and  d  -  73.4  ym. 
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Measured  times  for  70  and  90%  consumption  of  initial  agglomerate  mass. 
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largest  reduction  in  residence  time  is  observed  near  kinetic-controlled 
conditions.  This  behavior  agrees  with  the  trends  of  the  empirical 
parameters  as  du  is  varied  shown  in  Table  3,  i.e.,  the  a^  and  ?  in¬ 
crease  slightly  as  du  decreases,  tending  to  increase  rates  of  carbon 
consumption  at  kinetic-controlled  conditions.  In  general,  however,  these 
effects  are  small — roughly  comparable  to  present  experimental  accuracy. 

Effects  of  flame  conditions  and  initial  agglomerate  diameter  exert 
stronger  influences  on  residence  times  for  agglomerate  reaction  than  du- 
The  effect  of  flame  temperature  variations  for  <|>  *  0.6  and  1.0  is  illus¬ 
trated  in  Figs.  39  and  40.  Reacted  mass  fraction  is  plotted  as  a  function 

of  time  for  different  flame  temperatures  with  d  and  d  held  constant. 

po  u 

Predictions  for  the  models  using  both  reaction  mechanisms  are  shown  along 
with  the  measurements.  Increasing  temperature  results  in  appreciable  re¬ 
ductions  in  the  residence  time  required  to  react  a  given  fraction  of  the 
agglomerate  mass — particularly  at  conditions  approaching  reaction-control, 
e.g.,  <J>  -  1  in  Fig.  40.  Near  diffusion-controlled  conditions,  e.g.,  the 
high  temperature  results  for  <J>  ■  0.6  in  Fig.  39,  the  effect  of  temperature 
is  less  significant.  Here,  reaction  rates  are  primarily  dependent  on 
gas-phase  transport  properties  which  are  only  mildly  Influenced  by  tem¬ 
perature  changes  for  these  conditions. 

Similar  results  are  plotted  in  Figs.  41-43.  In  this  case,  <f>  is  varied 
for  a  constant  temperature  in  each  figure.  Over  this  test  range,  reducing 
4>  from  unity  results  in  significant  reduction  in  particle  lifetimes.  This 
occurs  since  oxygen  is  the  most  effective  reactant  of  carbon.  At  high 
equivalence  ratios,  there  is  little  oxygen  available  in  the  particle  sur¬ 
roundings  and  reaction  primarily  proceeds  via  OH — if  the  recent  results  of 
Neoh  et  al.  [18]  are  correct — and  OH  is  never  present  in  high  concentrations. 
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Figure  39.  Effect  of  ambient  temperature  on  the  variation 
of  reacted  mass  fraction  with  time.  $  -  0.6, 
d _  -  73. A  ym,  d  ■  300  nm. 
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Figure  40.  Effect  of  ambient  temperature  on  the  variation 

of  reacted  mass  fraction  with  time.  $  -  1.0, 

d  “73.4  ym,  d  ■  300  run. 
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Figure  41.  Effect  of  ambient  equivalence  ratio  on  the  variation 

of  reacted  mass  fraction  with  time.  T.  *  1665  K, 

d  -73.4  ym,  d  ■  300  nm. 
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Figure  42.  Effect  of  ambient  equivalence  ratio  on  the  variation 

of  reacted  mass  fraction  with  time.  ■  1815  K, 

d  *73. 4  um, d  ■  300  nm.  1 
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Figure  43.  Effect  of  ambient  equivalence  rati.'  un  ■  >*  • 
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The  effect  of  initial  agglomerate  size  on  residence  times  required  for 

carbon  reaction  is  illustrated  in  Fig.  44 ,  for  <f>  “  1  and  flame  temperatures 

of  1683  and  1952  K.  As  noted  earlier,  residence  times  for  reaction  are 

X  to  2 

proportional  to  for  the  low  Reynolds  numbers  considered  during 

present  tests.  Residence  time  requirements  are  nearly  proportional  to 
dp0  as  the  process  approaches  reaction-control.  This  condition  is  repre¬ 
sented  by  the  test  results  at  1683  K  in  Fig.  44.  Therefore,  improved 
atomization,  yielding  smaller  initial  agglomerate  diameters,  is  one  of 
the  most  effective  methods  for  improving  combustion  efficiency. 

4.2  Agglomerates  of  Carbon-Black  Blends 
4.2.1  Appearance 

Blends  of  carbon-blacks  having  different  ultimate  carbon 
particle  sizes  are  likely  for  formulations  of  practical  slurry  fuels,  in 
order  to  obtain  acceptable  rheological  and  stability  properties.  The  final 
phase  of  the  present  investigation  involved  a  brief  examination  of  the  com¬ 
bustion  properties  of  blends.  Work  was  limited  to  testing  and  analysis 
of  a  single  bimodal  blend  consisting  of  carbon-blacks  having  ultimate 
carbon  particle  sizes  of  70  and  300  nm,  with  equal  mass  loading  of  each 
size. 

The  appearance  of  agglomerates  from  the  bimodal  blend  was  similar  to 
the  monodlsperse  carbon-black  formulations  considered  during  this  study. 

The  main  difference  was  at  the  start  of  reaction,  where  the  bimodal 
agglomerates  were  somewhat  smoother  than  the  monodlsperse  agglomerates. 

This  follows  from  the  capability  of  the  smaller  particles  to  fit  in  the 
space  between  the  large  particles  in  the  blend. 


Figure  44.  Effect  of  initial  particle  diameter  on  the  variation 

of  reacted  mass  fraction  with  time.  $  ■  1,  d„  11  300  nm 


4.2.2  Empirical  Parameters 

Data  obtained  for  the  blmodal  blends  are  summarized  in 

Appendix  A.  The  initial  density  of  blmodal  agglomerates  was  roughly  16% 

3 

higher  than  the  monodisperse  agglomerates — 2108  kg/m  as  opposed  to 
3 

1820  kg/m  .  This  is  expected  from  solid-state  theory,  since  the  smaller 
particles  help  to  fill  in  the  voids  between  the  larger  particles. 

The  variation  of  blmodal  agglomerate  density  as  a  function  of  extent 
of  reaction  is  illustrated  in  Fig.  45  for  the  three  flame  conditions 
tested.  The  empirical  fit  developed  for  monodisperse  agglomerates, 

Eq.  (1),  is  also  shown  on  the  figure  for  comparison.  The  variation  in 
density  with  extent  of  reaction  is  relatively  independent  of  flame 
conditions.  The  density  of  the  blmodal  agglomerate  varied  in  the  same 
manner  as  for  the  monodisperse  agglomerates,  aside  from  the  different 
initial  value  of  density.  This  yields  the  following  expression  for  the 
density  of  the  blmodal  agglomerate 

pp  -  2108  (1  -  e)0,6  (25) 

These  results  imply  similar  behavior  for  the  development  of  the  pore 
structure  for  both  monodisperse  and  blmodal  agglomerates.  The  pore 
layer  has  a  characteristic  thickness  less  than  25%  of  the  particle  radius 
for  e  <  75%,  but  penetrates  to  the  center  of  the  agglomerate  during  the 
latter  stages  of  reaction. 

The  pore  structure  of  the  blmodal  agglomerate  enhances  convective 
transport,  however,  the  effect  was  less  pronounced  than  for  monodisperse 
agglomerates.  The  transport-enhancement  factors  for  the  blmodal  agglomerate 
are  illustrated  in  Fig.  46,  along  with  the  correlations  found  for  present 
tssts  with  monodisperse  agglomerates  and  the  earlier  findings  for  large 
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monodisperse  agglomerates  [6,7].  The  blends  behave  qualitatively  the 
same  as  the  monodisperse  agglomerates.  The  value  of  £  was  nearly  unity 
at  the  start  of  the  process.  As  the  pores  develop,  however,  £  increases 
reaching  a  maximum  near  e  “  0.75.  The  reduction  of  £  in  the  latter 
stages  of  the  process  corresponds  to  the  region  when  the  pore  layer  has 
reached  the  center  of  the  agglomerate.  The  variation  of  £  with  e  is 
essentially  independent  of  flame  conditions. 

The  main  difference  in  £  between  bimodal  and  monodisperse  agglomerates 
is  that  maximum  values  of  £  are  50Z  lower  for  bimodal  agglomerates.  It 
appears  that  the  smaller  particles  in  the  bimodal  agglomerate  tend  to 
inhibit  the  development  of  pores.  Therefore,  the  gas  is  less  able  to 
percolate  through  the  particle,  i.e.,  the  small  particles  tend  to  clog  the 
void  space  between  the  larger  particles. 

Area/reactivity  multiplication  factors  for  the  Oj-CO^-H  ,0  reaction 
mechanism  are  illustrated  in  Fig.  47.  Results  for  bimodal  agglomerates  as 
well  as  correlations  for  the  monodisperse  agglomerates  are  shown  on  the 
figure.  The  are  relatively  independent  of  flame  conditions  and  tend  to 
Increase  somewhat  with  increasing  e.  The  values  of  the  a^  for  the  bimodal 
agglomerates,  however,  are  roughly  50Z  of  the  values  found  for  monodisperse 
agglomerates.  Consistent  with  the  lower  value  of  £,  the  presence  of  the 
small  particles  appears  to  impede  the  development  (or  extent)  of  pores,  so 
that  less  active  surface  area  is  available  for  reaction. 

The  a^  for  the  OH-Oj  reaction  mechanism  are  illustrated  in  Fig.  48  for 
the  bimodal  blends  along  with  the  values  found  for  monodisperse  agglomerates. 
As  before,  the  a^  are  independent  of  flame  conditions  and  tend  to  increase 
as  e  increases  and  the  pore  layer  develops.  The  a^  for  the  bimodal  blend  are 
also  roughly  half  the  values  for  monodisperse  agglomerates. 


REACTIVITY 


oxygen 


reacted  mass  , 


The  empirical  parameters  for  the  blmodal  blends  were  correlated  as 

functions  of  £  and  d  ,  since  effects  of  flame  conditions  were  small. 

po 

These  correlations  are  summarized  in  Table  5.  The  present  data  were 

limited  to  a  single  initial  agglomerate  diameter.  Therefore,  the  effect 

of  d  has  arbitrarily  been  assumed  to  be  the  same  as  for  the  mono- 
po 

disperse  agglomerates.  In  addition  to  the  correlations,  constant  average 
values  of  the  empirical  parameters  are  listed  in  Table  5 — representative 
of  conditions  where  £  >  0.2.  The  average  values  found  for  the  mono- 
disperse  agglomerates  are  also  summarized  in  the  table  for  comparison.  As 
noted  earlier,  empirical  parameters  for  the  blmodal  blend  are  roughly 
50%  as.  large  as  the  values  for  monodisperse  agglomerates. 

4.2.3  Life  Histories 

The  burning  rate  results  for  the  blmodal  agglomerates 
were  similar  to  the  findings  for  monodisperse  agglomerates.  Use  of  the 
empirical  parameters  summarized  in  Table  5  yielded  good  agreement  between 
predictions  and  measurements.  The  use  of  the  constant  average  values  of 
the  parameters  resulted  in  only  a  slight  loss  in  the  accuracy  of  the  pre¬ 
dictions. 

Particle-life  histories  for  the  blmodal  agglomerate  are  Illustrated 
in  Figs.  49-51.  The  variation  of  particle  mass  and  diameter,  normalized 
by  their  initial  values,  and  the  particle  temperature,  normalized  by  the 
flame  temperature,  are  plotted  as  a  function  of  time  for  dpQ  ■  73.5  pm. 
Predictions  of  these  quantities,  for  both  reaction  mechanisms,  are  also 
shown  on  the  figures. 

Similar  to  the  earlier  results,  the  comparison  between  predictions  and 
measurements  in  Figs.  49-51  is  reasonably  good.  Aside  from  predicting 
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Table  5.  Suanary  of  Agglomerate  Structure  Parameters  (Bimodal  Blend) 
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Figure  49.  Particle-life  history  for  tp  -  0.6,  Tf  -  1663  K 
and  d  -  73.5  Mm  for  the  bloods!  agglomerate. 


burning  rates  reasonably  well,  which  is  expected,  the  analysis  yields  good 
predictions  of  particle  heat-up  and  maximum  temperature  levels. 

The  reduced  values  of  £  and  the  a^  for  the  blmodal  agglomerates  result 
in  somewhat  longer  residence  times  for  reaction.  Times  required  to  react 
60  and  90 %  of  the  initial  particle  mass  are  summarized  in  Table  6  for 
monodisperse  agglomerates  having  d^  -  70  and  300  nm  and  the  blmodal  agglom¬ 
erate  employing  these  constituents.  The  greatest  increase  in  residence 
time  is  observed  for  the  near  lcinetic-control  condition,  where  tgQg  is 
nearly  50%  longer  for  the  blend  than  for  the  monodisperse  agglomerates. 

The  effect  of  blending  is  more  complex  near  diffusion-control  con¬ 
ditions.  The  value  of  t^^  is  50%  longer  for  the  blend,  but,  the  blend 
seems  to  catch-up  during  the  latter  stages  of  the  process  and  t^^  is  only 
10-20%  longer  for  the  blend.  Comparing  the  particle-life  histories  for 
the  cases  in  Table  6  provides  a  clue  for  this  behavior.  The  agglomerate 
diameter  decreases  more  slowly  for  the  blend,  so  that  the  particle  is 
larger  near  the  end  of  the  process.  This  compensates  for  the  lower  value 
of  £  by  providing  more  apparent  surface  area. 

These  results  Indicate  that  residence  time  requirements  are  somewhat 
insensitive  to  changes  in  £  and  the  a^ — except  for  effects  of  a^  near 
kinetic-controlled  conditions.  This  insensitivity  is  due  to  the  relative 
complexity  of  the  reaction  mechanism,  so  that  changes  in  any  one  parameter 
tend  to  be  compensated  by  other  aspects  of  the  process. 

While  the  residence  times  for  the  blmodal  blend  are  somewhat  longer 
than  the  monodisperse  agglomerates,  for  a  given  initial  agglomerate  size, 
it  should  also  be  recalled  that  the  initial  density  of  the  blended  agglom¬ 
erate  is  16%  greater.  This  implies  that  agglomerates  of  blends  will 


Table  6.  Effect  of  Blending  on  Agglomerate  Combustion  Times 


Predicted  value  in  parenthesis. 


generally  be  smaller  than  those  of  monodisperse  carbon  blacks,  for  a  given 
slurry  drop  mass.  Since  smaller  particles  have  shorter  residence  time 
requirements,  this  compensates  for  the  somewhat  longer  residence  time 
requirements  of  blends.  Furthermore,  if  use  of  blends  improves  atomization 
properties  of  the  slurries,  smaller  drop  sizes  could  easily  yield  shorter 
residence  time  requirements  for  blends. 

5.  Conclusions 

The  present  model  provided  a  good  correlation  of  agglomerate  reaction 
rates  and  life  histories  using  either  the  Neoh  et  al.  [18]  or  extended 
Libby  and  Blake  [17]  carbon  reaction  mechanisms.  Combining  present  and 
earlier  [6,7]  results,  the  test  range  includes:  initial  agglomerate 
diameters  of  10-1000  ira,  agglomerate  Reynolds  numbers  of  0-200,  flame 
equivalence  ratios  of  0.2-1. 4,  C/8  ratios  (mass  basis)  of  2.25-3.10  and 
flame  temperatures  of  779-1953  K  at  atmospheric  pressure.  This  agreement 
is  encouraging,  removing  earlier  uncertainties  concerning  effects  of 
particle  size  as  well  as  flame  temperature  and  C/H  ratio  [6,7].  Furthermore, 
present  findings  support  the  conceptually  simple  carbon  reaction  mechanism 
of  Neoh  et  al.  [18]  at  high  fuel  equivalence  ratios. 

A  shrinking  sphere  model  was  used  in  the  analysis  which  Incorporates 
empirical  factors  to  treat  effects  of  pores  on  density,  enhancement  of 
convective  transport,  and  effective  surface  area  for  chemical  reactions. 

The  present  area/reactivity  multiplication  factors  can  also  Incorporate 
effects  of  catalyst  [6,7].  The  empirical  factors  were  relatively  inde¬ 
pendent  of  particle  size  and  flame  conditions.  Particle  density  decreased 
monotonically  with  increasing  extent  of  reaction.  Convection  enhancement 
was  Initially  relatively  small,  but  increased  as  the  pore  structure  de¬ 
veloped  which  allowed  the  flow  to  penetrate  the  porous  surface. 


Effects  of  ultimate  carbon  particle  size  of  the  carbon-black,  as 
well  as  a  limited  study  of  effects  of  blending  carbon  blacks  was  also 
considered.  Reducing  the  ultimate  carbon  particle  size  from  300  to  70  nm 
resulted  in  slightly  higher  burning  rates,  yielding  reductions  in  resi¬ 
dence  times  for  particle  reaction  on  the  order  of  10-20Z — depending  on 
flame  conditions. 

Only  a  bimodal  blend  of  carbon  black  having  ultimate  carbon  particle 
sizes  of  70  and  300  nm  (50%  by  weight  of  each)  was  considered.  This  blend 
had  an  initial  agglomerate  density  that  was  16%  higher  than  the  mono- 
disperse  agglomerates.  Residence  times  needed  to  react  given  fractions 
of  the  carbon  mass  for  the  blend  were  longer  than  for  monodisperse  agglom¬ 
erates  having  the  same  initial  diameter.  Residence  time  increases  ranged 
from  10-20%  at  diffusion-controlled  conditions  to  50%  at  kinetlc-controlled 
conditions.  The  increased  mass  of  carbon  to  be  reacted  in  the  blend  accounts 
for  a  portion  of  this  increase.  The  effect  of  smaller  particles  filling 
void  spaces  between  larger  particles,  thus  inhibiting  the  development  of 
pores,  is  probably  also  a  factor  contributing  to  increased  residence  time 
requirements  for  blends.  These  effects  also  caused  changes  in  the  empirical 
parameters  found  for  the  blends.  Transport  enhancement  and  area/reactivity 
multiplication  factors  for  blends  were  roughly  50%  of  the  values  found  for 
monodisperse  agglomerates,  indicating  some  suppression  of  pore  development 
for  blends. 

The  present  results  indicate  that  carbon-black  formulation  affects 
agglomerate  combustion  properties.  Use  of  carbon-blacks  having  smaller 
ultimate  carbon  particle  sizes  tends  to  reduce  residence  time  require¬ 
ments.  Blends  require  somewhat  longer  residence  times  for  reaction  than 
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monodisperse  agglomerates  of  chair  conaCltuants — for  cha  same  Initial 
agglomerate  size.  These  af facta  ware  not  large,  however,  particularly 
at  flame  conditions  where  agglomerate  reaction  rates  were  high.  There¬ 
fore,  if  blending  yields  improved  atomisation  properties,  the  production 
of  smaller  drops  could  easily  result  in  shorter  residence  time  require¬ 
ments  for  slurry  blends  during  combustor  operation. 

Present  findings  were  limited  to  a  single  biaodal  blend.  In  order 
to  establish  the  combustion  properties  of  blends,  other  ultimate  carbon 
particle  sizes  and  variations  in  the  proportions  of  each  size  should  be 
considered.  It  should  also  be  noted  that  residence  times  found  here  are 
for  fixed  ambient  conditions,  while  agglomerates  in  practical  combustors 
will  experience  a  range  of  conditions  during  their  lifetimes.  The  present 
analysis  provides  a  means  of  estimating  effects  of  variations  in  ambient 
conditions — at  least  over  the  range  of  conditions  where  it  has  been  evalu¬ 
ated.  In  particular,  evaluation  thus  far  has  been  limited  to  atmospheric 
pressure;  therefore,  capabilities  for  predicting  effects  of  pressure  are 


uncertain. 


Appendix  A 


Summary  of  Data 


A, 

A.l  Particle  Environment 


Gas 

Temperature  and 

Composition 

,  <j>  -  0.2 

X 

(mm) 

T 

<K) 

Mole  Fraction 

°2 

C02 

CO 

h2° 

«2 

H2 

OH 

10 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691 

— 

0.0004 

20 

1668 

0.5314 

0.0664 

— 

0.1328 

0.2691 

— 

0.0004 

30 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691 

— 

0.0004 

40 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691 

— 

0.0004 

50 

1667 

0.5314 

0.0664 

— 

0.1328 

0.2691 

— 

0.0004 

Gas 

Temperature  and  Composition 

,  $  -  0.6 

X 

(mm) 

T 

<K) 

Mole  Fraction 

°2 

co2 

CO 

V 

H2 

H2 

OH 

5 

1662 

0.0860 

0.0637 

— 

0.1293 

0.7208 

— 

0.0002 

10 

1663 

0.0858 

0.0640 

— 

0.1291 

0.7208 

— 

0.0003 

20 

1663 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

30 

1664 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

40 

1663 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

50 

1664 

0.0858 

0.0644 

— 

0.1291 

0.7208 

— 

0.0003 

A0H  mole  fraction  was  obtained  from  CEC-72  equilibrium  calculation. 
^0  mole  fraction  was  obtained  from  H  and  0  element  balances. 
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Gas  Temperature  and  Composition,  $  -  0.6 


X 

T 

Mole  Fraction 

(mm) 

(K) 

°2 

C°2 

CO 

V 

N2 

"2 

j 

10 

1810 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

r* 

20 

1814 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

30 

1815 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

** 

40 

1814 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

i 

50 

1814 

0.0973 

0.0730 

0.0001 

0.1458 

0.6832 

— 

x  T 
(nm)  (K) 


Gas  Temperature  and  Composition,  $  »  0.6 


_ Mole  Fraction _ 

C0„  CO  H„0  N. 


Gas  Temperature  and  Composition,  $  *  1.0 


x  T  _ Mole  Fraction 


(mm) 

<K) 

°2 

C02 

CO 

h2o 

N2 

H2 

OH 

5 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

10 

1683 

0.0002 

0.0645 

0.0003 

0.1294 

0.8054 

0.0002 

20 

1684 

0.0002 

0.0645 

0.0003 

0.1294 

0.8054 

0.0002 

0.0001 

30 

1683 

0.0002 

0  0645 

0.0003 

0.1295 

0.8054 

0.0002 

40 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

50 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

0.0001 

60 

1683 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

70 

1681 

0.0002 

0.0645 

0.0003 

0.1295 

0.8054 

0.0002 

80 

1678 

0.0002 

0.0646 

0.0003 

0.1295 

0.8054 

0.0002 

90 

1675 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0002 

0.0001 

100 

1668 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0002 

0.0001 

110 

1665 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0002 

0.0001 

120 

1661 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0001 

0.0001 

130 

1657 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0001 

0.0001 

140 

1648 

0.0002 

0.0646 

0.0002 

0.1295 

0.8054 

0.0001 

0.0001 

150 

1640 

0.0001 

0.0646 

0.0002 

0.1295 

0.8055 

0.0001 

160 

1625 

0.0001 

0.0646 

0.0002 

0.1295 

0.8055 

0.0001 

0.0 

170 

1606 

0.0001 

0.0647 

0.0001 

0.1295 

0.8055 

0.0001 

0.0 

180 

1597 

0.0001 

0.0647 

0.0001 

0.1295 

0.8055 

0.0001 

0.0 

190 

1573 

0.0001 

0.0647 

0.0001 

0.1296 

0.8055 

0.0001 

0.0 

200 

1564 

0.0001 

0.0647 

0.0001 

0.1296 

0.8055 

0.0001 

0.0 

Gas  Temperature  and  Composition,  <J>  •  1.0 


10 

1816 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

20 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

30 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

40 

1820 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

50 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

60 

1819 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

70 

1817 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

30 

1817 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

90 

1815 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

100 

1812 

0.0005 

0.0722 

0.0008 

0.1455 

0.7805 

0.0004 

no 

1809 

0.0005 

0.0722 

0.0007 

0.1456 

0.7806 

0.0004 

120 

1805 

0.0005 

0.0723 

0.0007 

0.1455 

0.7806 

0.0004 

130 

1801 

0.0005 

0.0723 

0.0007 

0.1455 

0.7806 

0.0004 

140 

1793 

0.0005 

0.0723 

0.0006 

0.1456 

0.7806 

0.0003 

150 

1785 

0.0009 

0.0724 

0.0006 

0.1456 

0.7806 

0.0003 

160 

1773 

0.0004 

0.0724 

0.0006 

0.1457 

0.7806 

0.0003 

170 

1754 

0.0004 

0.0725 

0.0005 

0.1456 

0.7807 

0.0003 

180 

1742 

0.0003 

0.0725 

0.0004 

0.1458 

0.7807 

0.0003 

190 

1730 

0.0003 

0.0726 

0.0004 

0.1457 

0.7808 

0.0002 

200 

1718 

0.0003 

0.0726 

0.0004 

0.1457 

0.7808 

0.0002 

0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 
0.0001 
0.0001 
0.0001 
0.0001 
0.0001 


[• 


Gas  Temperature  and  Composition,  <J>  ■  1.0 


x  T 


Mole  Fraction 


(mm) 

(K) 

°2 

n 

o 

N> 

CO 

V 

N2 

H2 

10 

1949 

0.0013 

0.0806 

0.0019 

0.1641 

0.7509 

0.0009 

20 

1952 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

0.0009 

30 

1953 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

0.0009 

40 

1952 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

0.0009 

50 

1952 

0.0013 

0.0806 

0.0020 

0.1641 

0.7509 

0.0009 

60 

1951 

0.0013 

0.0806 

0.0019 

0.1642 

0.7509 

0.0009 

70 

1950 

0.0013 

0.0806 

0.0019 

0.1642 

0.7509 

0.0009 

30 

1948 

0.0013 

0.0806 

0.0019 

0.1642 

0.7509 

0.0009 

90 

1946 

0.0012 

0.0806 

0.0019 

0.1644 

0.7509 

0.0009 

100 

1944 

0.0012 

0.0807 

0.0019 

0.1643 

0.7509 

0.0009 

110 

1940 

0.0012 

0.0807 

0.0018 

0.1643 

0.7510 

0.0009 

120 

1935 

0.0012 

0.0808 

0.0018 

0.1643 

0.7510 

0.0008 

130 

1930 

0.0011 

0.0808 

0.0017 

0.1644 

0.7510 

0.0009 

140 

1924 

0.0011 

0.0809 

0.0016 

0.1644 

0.7511 

0.0008 

150 

1913 

0.0010 

0.0810 

0.0015 

0.1745 

0.7511 

0.0007 

160 

1901 

0.0010 

0.0811 

0.0014 

0.1645 

0.7512 

0.0007 

170 

1891 

0.0009 

0.0812 

0.0013 

0.1645 

0.7512 

0.0007 

130 

1879 

0.0008 

0.0813 

0.0012 

0.1647 

0.0006 

0.0006 

190 

1860 

0.0008 

0.0814 

0.0011 

0.1647 

0.7514 

0.0006 

200 

1851 

0.0007 

0.815 

0.0010 

0.1648 

0.7614 

0.0006 

OH 


0.0003 

0.0003 
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Gas  Temperature  and  Composition,  $  •  1.4 


x  T _ Mole  Fraction 


(mm) 

(K) 

0o 

a* 

C02 

CO 

V 

*2 

H2 

OH 

5 

1622 

— 

0.0421 

0.0534 

0.1432 

0.7092 

0.0542 

— 

10 

1624 

— 

0.0427 

0.0540 

0.1412 

0.7092 

0.0563 
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1658 

1.523 

0.134 

182. 7 

60.7 

130.7 

1655 

1.548 

0.145 

59.8 

119.8 

1.550 

0.158 

Particle-Life  History,  $ 


1663  K 


X 

d 

m 

T 

u 

t 

P 

P 

P 

P 

(mm) 

(pm) 

(ng) 

(K) 

(m/s) 

(8) 

0.0 

73.5 

397.7 

607 

0.435 

0.0 

2.0 

71.7 

331.7 

1998 

0.581 

0.0038 

3.2 

68.6 

235.8 

2033 

0.628 

0.0058 

4.4 

65.4 

166.1 

2044 

0.687 

0.0077 

7.6 

59.2 

102.0 

2057 

0.721 

0.012 

10.5 

54.5 

60.3 

2069 

0.732 

0.016 

12.8 

49.2 

29.7 

2084 

0.763 

0.019 

16.7 

43.4 

13.9 

2092 

0.789 

0.024 

19.1 

39.4 

8.4 

2101 

0.802 

0.027 

23.2 

32.6 

5.6 

2109 

0.837 

0.032 

26.7 

36.7 

2.5 

2118 

0.871 

0.036 

30.2 

19.1 

0.7 

2138 

0.888 

0.040 

33.9 

12.3 

0.006 

2126 

0.912 

0.044 

37.6 

6.9 

0.0 

2098 

0.997 

0.048 

39.7 

2.4 

0.0 

2079 

1.073 

0.050 

0.0 

73.5 

404.3 

610 

1.5 

68.7 

246.8 

2174 

2.4 

64.6 

155.4 

2191 

3.0 

62.5 

121.4 

2198 

4.1 

59.5 

75.5 

2211 

5.2 

55.6 

51.9 

2223 

6.2 

49.7 

31.3 

2229 

7.3 

48.1 

23.6 

2236 

8.5 

41.5 

9.43 

2242 

0.1 

35.2 

4.42 

2257 

.1.3 

30.8 

3.14 

2259 

.2.4 

26.6 

2.47 

2261 

.4.3 

19.6 

0.52 

2262 

.6.7 

1.03 

0.004 

2262 

.8.0 

5.7 

0.0 

2260 

.9.6 

1.8 

0.0 

2158 

0.0037 

0.0050 

0.0057 

0.0071 

0.0085 

0.0098 

0.0115 

0.0128 

0.0145 

0.0158 

0.0170 

0.0187 

0.0209 

0.0221 

0.0234 
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Psrticle-Lifs  History,  <j>  -  1.0,  Tf  -  1683  K 


X 

d 

a 

T 

u 

t 

P 

P 

P 

P 

(am) 

(pm) 

(ng) 

<K) 

(«/•) 

(a) 

0.0 

73.4 

402.1 

609 

0.441 

0.0 

1.9 

73.4 

401.0 

1516 

0.583 

0.0037 

2.8 

73.3 

397.4 

1605 

0.617 

0.0052 

4.2 

73.2 

388.9 

1640 

0.664 

0.0074 

6.1 

72.9 

377.0 

1647 

0.720 

0.0101 

8.2 

72.5 

364.5 

1648 

0.775 

0.0129 

15.8 

71.6 

328.8 

1648 

0.929 

0.0218 

21.3 

70.9 

307.5 

1650 

1.015 

0.0275 

30.5 

69.8 

278.1 

1650 

1.126 

0.0360 

41.8 

68.9 

248.3 

1650 

1.228 

0.0456 

53.0 

68.0 

223.8 

1651 

1.302 

0.0544 

64.2 

67.1 

202.4 

1651 

1.361 

0.0628 

68.9 

66.7 

194.2 

1652 

1.382 

0.0662 

79.8 

65.2 

184.9 

1654 

1.393 

0.0741 

91.9 

63.9 

168.8 

1655 

1.402 

0.0827 

105.5 

63.1 

160.6 

1654 

1.427 

0.0924 

122.5 

61.6 

140.7 

1651 

1.473 

0.104 

158.0 

60.3 

130.4 

1650 

1.493 

0.128 

171.5 

57.4 

107.6 

1647 

1.510 

0.137 

192.8 

56.1 

92.5 

1641 

1.539 

0.151 

215.2 

55.8 

81.3 

1632 

1.547 

0.168 
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Particle-Life  History,  <p  m  1.0,  -  1683  K 


X 

d 

m 

T 

u 

t 

P 

P 

P 

P 

(mm) 

<V») 

(ng) 

(K) 

(m/s) 

(s) 

m m 

73.3 

434.7 

609 

0.442 

0.0 

1.8 

73.3 

434.5 

1441 

0.565 

0.0035 

2.9 

73.3 

433.3 

1584 

0.604 

0.0053 

4.2 

73.2 

430.3 

1643 

0.641 

0.0074 

6.3 

73.1 

424.2 

1664 

0.695 

0.0106 

7.9 

73.0 

419.5 

1667 

0.731 

0.0129 

15.5 

72.5 

400.5 

1668 

0.863 

0.022 

22.6 

72.1 

385.4 

1668 

0.954 

0.030 

31.9 

71.7 

368.3 

1667 

1.045 

0.039 

41.4 

71.3 

352.9 

1668 

1.117 

0.048 

52.7 

70.8 

336.5 

1668 

1.181 

0.058 

63.2 

70.4 

322.4 

1668 

1.229 

0.067 

73.0 

70.0 

310.5 

1669 

1.265 

0.075 

91.4 

69.4 

289.7 

1668 

1.317 

0.089 

104.2 

69.0 

276.6 

1669 

1.345 

0.098 

139.9 

66.6 

252.1 

1668 

1.402 

0.124 

158.2 

65.7 

233.9 

1661 

1.409 

0.137 

169.4 

63.5 

213.0 

1657 

1.419 

0.145 

187.8 

62.7 

201.7 

1649 

1.431 

0.158 

207.6 

60.8 

182.6 

1633 

1.454 

0.173 
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